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ABSTRACT 

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active 
galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near 
Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP-i-GOODS-Herschel survey. We derive 
detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we 
then compare to those of the AGNs. At z ~ 1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as 
well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In 
contrast, at z ~ 2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, 
which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests 
of both these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution 
among AGNs and inactive galaxies are similar at z > 1.5. At z < 1, however, clear and significant enhancements are seen in the 
SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and 

host nrnnertiYs with reHshift tnwarrls; pi minor rolo for moraors; pmH intorpirtions pit r > 1 5 


1. Introduction 

The study of galaxy structure (or morphology, used interchange- 
ably in this work) has provided important insights into galaxy 
formation and evolution and is a key observable in constraints 
on theoretical models. It is conclusively linked to other important 
galaxy properties such as mass, baryonic content, star-formation 
history, interaction state and environment (e.g.. Dressier 1980; 
Roberts & Haynes 1994; Kennicutt 1998; Strateva et al. 2001; 
Wuyts et al. 2011). In this sense, structure is a sensitive mea- 
sure of the history of galaxy growth and can be used to constrain 
galaxy evolution models. Another important tracer of the evolu- 
tionary state of a galaxy is the star-formation rate (SFR). This 
is sensitive to the gas content and infall onto a galaxy, as well 
as the processes that inject energy and momentum into the ISM 
and regulate the overall efficiency of cold gas fragmentation and 
the formation of new stars. 

Active Galactic Nuclei (AGNs) are found in a subset of 
galaxies where a central supermassive black hole (SMBH) is 
accreting material at a sufficient level to be detectable us- 
ing many characteristic tracers of high-energy activity, such as 
strong X-ray, non-thermal radio, hot dust or high excitation line 
emission. There are many theories as to what triggers and fu- 
els SMBH accretion (e.g., Hopkins et al. 2005; Fanidakis et al. 
2012; Hirschmann et al. 2012; see Alexander & Hickox 2012 
for a contemporary discussion). These may be broadly divided 
into secular processes (bars and spirals, possibly minor mergers) 
that take several galaxy dynamical timescales {tdyn) to bring gas 
to the nucleus from galaxy scales, or “violent” processes (ha- 
rassment, major mergers, violent disk instability), which operate 
through varying, typically external, torques that change on tdyn, 
and are expected to circumvent secular inflow. Since tdyn is the 
characteristic duration over which a galaxy responds to a violent 
disturbance and relaxes into a new state of dynamical equilib- 
rium, the relationship between galaxy structure and the occur- 
rence or strength of nuclear activity is a vital indicator of the 
relative importance of violent and secular AGN fueling mecha- 
nisms. 

Star-formation is also affected differently by these two 
classes of processes. Stable galactic gas disks form stars with 
a relatively low efficiency, which current studies suggest is 
modulated by the interplay of dense gas and stellar feedback 
(Kim et al. 2013; Hopkins et al. 2013a). Violent processes lead 
to the strong infiow of gas into the centers of galactic po- 
tential wells (Sanders et al. 1988) or into dense small-scale 
clumps (Dekel et al. 2009), both of which result in star- forming 
environments with an enhanced efficiency of star-formation 
(Genzeletal. 2010; Daddietal. 2010b). Observationally, star- 
forming galaxies show a clear correlation between SFR and stel- 
lar mass (M,) (Noeske et al. 2007; Elbaz et al. 2007; Daddi et al. 
2007; Santinietal. 2009; Rodighiero et al. 2011), sometimes 
called the star-formation (SF) sequence. The existence of this 
sequence implies that secular processes govern most star- 
forming galaxies and, indeed, structural studies of galaxies 
across much of cosmic time show that galaxies on the SF se- 
quence are disk-dominated, indicative of low levels of dynamical 
disturbance (Forster Schreiber et al. 2009; Wuyts et al. 2011). 
Starbursts are defined to have an elevated SFR, placing them 
above the SF sequence (Rodighiero et al. 2011; Sargent etal. 
2012), and are structurally consistent with being recently dis- 
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turbed (Wuyts et al. 2011; Kartaltepe et al. 2012). This associa- 
tion between structure and SF patterns highlights alternate path- 
ways for violent and secular processes, in terms of the connec- 
tion between SF and the morphological change in galaxies. 

The morphologies of AGN hosts have been extensively 
studied in various different ways, across a wide range of 
redshifts. In general, these studies have either used analytic 
measures of structure, such as light profile fitting, pixel dis- 
tribution statistics or moment measures (Grogin et al. 2003, 
2005; Pierce et al. 2007; Gabor et al. 2009; Simmons et al. 
2011; Schawinski et al. 2011; Bohm et al. 2013; Villforth et al. 
2014) or visual classification using a number of special- 
purpose schemes (Schawinski et al. 2007, 2010; Cistemas et al. 
2011; Kocevski et al. 2012). Some studies of more luminous 
optically-unobscured AGNs have also been undertaken, involv- 
ing the careful subtraction of nuclear point sources to re- 
veal the host galaxy more clearly (e.g., Bahcall et al. 1997; 
McEeod & McEeod 2001; Dunlop et al. 2003; Jahnke et al. 
2004; Sanchez et al. 2004; Guyon et al. 2006; Veilleux et al. 
2009). AGN hosts span a range of morphologies and levels of 
galaxy disturbance, with a sizeable fraction found in massive 
galaxies with substantial disks and low levels of merger activity. 
The most luminous systems tend to be in more bulge-dominated 
or early-type hosts. 

In parallel, the SF properties of AGN hosts is topic of 
much current work. Various tracers of the SFR have been ap- 
plied to AGNs, though traditional ones such as emission lines, 
UV continuum or optical spectral features can be strongly con- 
taminated by nuclear light, especially among luminous Type I 
AGNs. Despite the limited sensitivity of current datasets, far- 
infrared (FIR) wavelengths (> 50 pm) offer the best discrim- 
ination between emission from star-formation heated dust and 
AGN-heated dust, since the latter almost always exhibits a 
warmer distribution of temperatures, peaking in luminosity at 
~ 20 pm with a sharp drop off to the FIR (Netzer et al. 2007; 
Mullaney et al. 2011; Rosario et al. 2012). Herschel-based SFRs 
of active galaxies across a range of luminosities have been shown 
to agree reasonably well with those of massive inactive (i.e., non- 
AGN) star-forming galaxies (Mullaney et al. 2012; Santini et al. 
2012; Rovilos et al. 2012; Rosario et al. 2013b,c), but there is a 
strong tendency for X-ray bright AGNs to be found in SF hosts 
(Rosario et al. 2013b). Therefore, AGN hosts are primarily mas- 
sive galaxies forming stars at a normal rate given their stellar 
content. 

In this work, we present an exploration into connections be- 
tween the SFR and structure of AGN hosts at intermediate to 
high redshifts (z = [0.5, 2.5]), towards constraining the degree 
to which AGNs exhibit the relationships known among inactive 
galaxies (i.e., those without substantial on-going SMBH accre- 
tion). The primary goal of this paper is to lay the groundwork 
for a more extensive and detailed analysis of such relationships 
using forthcoming datasets spanning a number of important ex- 
tragalactic fields. Despite this, our work is already unique in its 
approach of testing SFR-structure trends in systems with signif- 
icant AGN activity across an important range in redshifts. Our 
study builds on the best available deep extragalactic data at X- 
ray, optical, NIR and FIR wavelengths, and leverages the sub- 
stantial efforts of the extragalactic community towards building 
coherent datasets in the GOODS fields. Datasets, selections and 
methods are introduced in Section 2. We discuss specific struc- 
tural patterns relevant to the galaxy population at these redshifts 
and compare them to the AGNs, revealing any differences be- 
tween these two populations (Section 3). In Section 4, we inves- 
tigate trends between SFR and structure, again comparing the 
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AGNs to the galaxy population. Finally, in Section 5, we discuss 
our results in the context of earlier work and weigh them quali- 
tatively against the predictions of key AGN triggering scenarios. 

We adopt a standard A-CDM Concordance cosmology, with 
Hq - 10 km s"^ Mpc“^ and Qa = 0.7. Stellar masses in this 
study are estimated assuming a Salpeter Initial Mass Function. 


2. Sample Selection, Datasets and Methods 

In this work, we concentrate on X-ray selected AGN, since lu- 
minous X-ray emission is an unambiguous hallmark of SMBH 
accretion activity. However, studies of the AGN population and 
the X-ray background suggest that X-ray selection systemati- 
cally misses the most obscured AGNs. Our conclusions from this 
work are strictly applicable only to such AGNs as are traced by 
X-ray selection. This caveat must be borne in mind when setting 
this work into the overall context of AGN demographics. 

2. 7. CAN DELS HST imaging 

A panchromatic spatially resolved view of substantial numbers 
of galaxies in the distant Universe is now available through 
the Cosmic Near-Infrared Deep Extragalactic Legacy Survey 
(CANDELS), which consists of HST ACS and WFC3 imag- 
ing of five separate deep extragalactic fields in multiple broad 
photometric bands spanning the observed UV to the near-IR. 
Full details of the design, data reduction and available for- 
mats of the survey are discussed in Groginetal. (2011) and 
Koekemoer et al. (2011). The great leap forward for galaxy 
structure work made possible by CANDELS comes from deep 
near- infrared imaging from the IR channel of the WFC3 camera. 
The reddest CANDELS band is F160W, which approximately 
corresponds to the H-band at 1.6 /urn. For galaxies at z = 2.5, 
F160W lies completely redward of rest- frame 4000 A. At these 
wavelengths, the continuum light from galaxies is representative 
of the total stellar content of the systems, and is significantly 
less biased by rest- frame ultra-violet light from young stars that 
can completely dominate in the observed optical and UV wave- 
lengths. Earlier studies of galaxy structure over large galaxy 
samples at HST resolution were restricted to observed wave- 
length < 1 /urn. With its capacity to capture the true stellar distri- 
bution of statistically complete samples of galaxies, CANDELS 
has greatly expanded the abilities of galaxy structural analysis at 
z ~ 2. 

We measure structure from the entire CANDELS mo- 
saics in the two fields associated with the Great Observatories 
Origins Deep Survey (GOODS Giavalisco et al. 2004). The 
CANDELS mosaic in GOODS-S consists of an ultradeep re- 
gion (Hubble Ultra-Deep field; 2' x 2.'3) enclosed by a deep re- 
gion (CANDELS-DEEP; 6.'8 x 10'), both of which are enclosed 
by a larger shallower region (CANDELS-Wide; 10' x 16'). The 
CANDELS mosaic in GOODS-N is similar in size to that in 
GOODS-S, with Deep and Wide subregions but without an ul- 
tradeep section. 

The HST bands cover different rest-frame wavelengths 
across the redshifts considered in this work. Therefore, it is im- 
portant to keep in mind that that when we discuss or highlight 
evolution in structural properties obtained from our own analy- 
sis, we specifically mean apparent evolution which may be sub- 
ject to morphological K-corrections. The main conclusions of 
the paper are based on structural and SFR comparisons between 
galaxy populations at similar redshifts, which are not subject to 
K-corrections. As a rough gauge for the reader, El 60 W (H-band) 


covers a rest-frame band at ~ SOOOA at z = 1 and ~ 5200A at 
z = 2. Two other HST optical bands relevant for appreciating the 
following analysis are F606W (V-band; ~ 3 000 A at z = 1 and 
^ 2000 A at z = 2) and F850LP (z-band; ~ 4500A at z = 1 and 
~ 3000Aatz = 2). 

2.1.1. GALFIT models of galaxy structure 

For our primary measure of galaxy structure in the two fields, 
we rely on light profile fits to the F160W (H-band) images 
of galaxies in CANDELS, performed using the GALFIT code 
(Pengetal. 2010). Details of the setup and fitting procedure, 
including the segmentation of the CANDELS mosaics, object 
identification, data preparations and initializations, can be found 
in van der Wei et al. (2012). 

Galaxies are fit using elliptical Sersic light profile models, 
yielding an estimate of the best-fit major-axis half-light radius 
Re, the Sersic index (n) and the ellipse axis ratio (q) for each 
galaxy through GALFIT. The Sersic index is the principal pa- 
rameter that governs the normalised light profile of galaxy, and, 
in a broad sense, is related to the dynamical state of the stellar 
matter (though the exact nature of these relationships and their 
evolution with redshift is yet unclear). The Sersic index is per- 
mitted to vary in the range of 0.5 (sub-exponential) to 8.0 (super- 
De Vaucouleurs). To ensure the fidelity of the fits, we restrict our 
study to sources with a F160W magnitude mh < 24.0 and with 
a GALFIT fiag= 0, indicating a sensible fit. 

We also include a set of GALFIT fits to CANDELS F160W 
and GOODS-S v2.0 F850LP (z-band) images of a subset 
of objects in the CANDELS GOODS-S field, those which 
overlap with the FIREWORKS multi-wavelength compilation 
(Wuyts et al. 2008). The GALFIT setup used for these fits is 
documented in Wuyts et al. (2011) and is very similar to that 
used for the H-band fits. We compared fits made using both se- 
tups for the FIREWORKS subset and found good consistency, 
with a median difference in n of 0.04 and a scatter in n of 0.25. 
About 15% differ by Aw = 1, typically towards a higher n using 
the setup from van der Wei et al. (2012) and mostly among high 
Sersic galaxies which have larger uncertainties in n. 

For this work, we fit galaxies, whether active or inactive, us- 
ing pure galaxy light profile models, without any special compo- 
nents to account for possible AGN contamination. Studies have 
shown that excess nuclear light from the AGN in the rest- frame 
optical bands at the level of tens of percent can systematically 
alter single profile galaxy fits, as well as other analytic mea- 
sures of structure. Therefore, earlier studies have tended to in- 
clude additional point source components to AGN profile fits, 
on the principle that such fits can distinguish between emission 
from a nuclear point source and the light of central bulges (which 
may appear very similar). There has been some calibration of 
this technique for AGN hosts imaged with the HST/ACS cam- 
era (Simmons & Urry 2008; Gabor et al. 2009). Even with the 
narrow ACS PSF, such studies found that this method, while 
generally sound, overestimates the nuclear point source frac- 
tion in a fraction of early-type systems, leading to an enhanced 
AGN component at the expense of a residual galaxy component 
with a lowered Sersic index. With the broader and more complex 
HST/WFC3-IR PSF, it is expected that these systematics would 
be more pronounced. To avoid being affected by fitting based 
biases, we employ only single component fits for all galaxies 
and use various tests to check for the effects of AGN contamina- 
tion when presenting results based on the GALFIT fits (Section 
3.3.1). 
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Scotter < 0,25 ■ 0,25 < Scotter < 0,50 ■ 0,50 < Scotter < 0,75 ■ Scotter > 0,75 



0,0 l _4 

0,00 0,25 0,50 0,75 1,00 0,00 0,25 0,50 0,75 1,00 

Final Interaction Metric Final Clumpiness Metric 

Fig. 1. The breakup of the mean seatter of visual elassifieation metries for galaxies at 0.5 < z < 2.5 binned by their final mean metrie 
value. The mean seatter for a galaxy is defined as the average absolute value of the differenee of the metrie between pairs of elassifiers 
for that galaxy. Galaxies with larger seatter have a larger level of disagreement between the individual visual elassifieations. The 
relative number of galaxies in four ranges of seatter are shown, represented by lighter to darker grayseale towards inereasing seatter. 
Galaxies are binned by the final interaetion metrie (IM, left panel) and the final elumpiness metrie (CM, right panel). 


2.1.2. Visual classification of galaxy structure and 
disturbance 

In the CANDELS GOODS-S field, all galaxies with mh < 
24.5 were examined by at least three human elassifiers, as part 
of a large on-going program of visual elassifieation by the 
CANDELS team. Eaeh elassifier looked at four images of the 
galaxy - 2-orbit F160W and F125W images, as well as F606W 
and F850LP images from the GOODS-S v2.0 release - along 
with its SExtraetor segmentation map on the F160W mosaie. 
Several details regarding a galaxy’s visual shape, eolor, asym- 
metry, elumpiness and disturbanee were noted in a systematie 
fashion using a GUI tool. A set of metries were used to order 
the various elassifieations down to a final redueed set, taking 
into aeeount the averaged reliability of elassifiers. More details 
about the plethora of visual elassifieation outputs, the methodol- 
ogy of elassifieation and a diseussion of elassifier reliability ean 
be found in Kartaltepe et al. (2014). 

In this work, we rely on visual estimates of the degree of 
gross galaxy disturbanee, using two different benehmarks. One 
is the “interaetion metrie” {IM), a general measure of the level 
of interaetion or merging in a galaxy. The CANDELS visual 
elassifieation tool gives elassifiers five ehoiees for eaeh exam- 
ined objeet, along a sequenee of inereasing apparent degree of 
interaetion. We assign a numerieal value to eaeh ehoiee, whieh 
enables us to easily eombine the outputs of different elassifiers. 
Objeets whieh are elearly undisturbed, with no obvious nearby 
neighbours and no signature of any interaetion, are given an 
IM - 0. At the other extreme, objeets with IM - 1 are in ob- 
vious late-stage mergers, with highly disturbed strueture, strong 
asymmetries, mueh elumpiness, frequently showing tidal tails 
and/or double or multiple nuelei. IM = 0.25 denote objeets in 
an apparent pair or multiple system separated by a few to sev- 


eral areseeonds with no elear signs of interaetion; these may be 
assoeiated but may also be line-of-sight alignments. IM =0.5 
are given to objeets with visual signatures of interaetion with 
other galaxies outside their H-band segmentation maps - these 
may be pre-mergers or elose eneounters. IM = 0.75 are objeets 
that have another galaxy within their segmentation areas and are 
probably in an interaeting state, but still show distinet strueture, 
likely beeause the merger is still at an early stage. The deeisions 
of individual reliable elassifiers may differ, but we average their 
IM values leading to final /Ms that eould be intermediate to the 
five prineipal assignments deseribed above. The visual IM for 
distant galaxies is most sensitive to major mergers, i.e., those 
where the optieal luminosity ratio of the eomponents ~ 1, or, 
alternatively, where the late-stage merger shows strong signs of 
disturbanee. Minor mergers will only be deteetable among the 
few big bright relatively loeal galaxies in the CANDELS fields, 
none of whieh are eonsidered in this work. 

To develop a sense of the level of disagreement in the visual 
assessment of interaetion aeross different elassifiers, we eonsider 
the seatter of IM for all galaxies from CANDELS/GOODS-S 
in the redshift interval 0.5 < z < 2.5. For eaeh galaxy with N 
elassifieations, we ealeulate the mean absolute value of the dif- 
ferenee of IM between the N{N - l)/2 pairs of elassifieations. 
Sinee the number of elassifiers for a galaxy eould be as low as 
= 4, this quantity is not neeessarily a statistieally aeeurate 
uneertainty on the final metrie, but eonsidering a large number 
of galaxies together, its distribution allows a simple measure of 
the variation in IM. We show this distribution in the left panel 
of Figure 1, splitting galaxies into 4 bins of their final eombined 
IM and eonsidering four bins in the mean seatter. It is immedi- 
ately apparent that galaxies elassified as undisturbed or strongly 
interaeting have only a small seatter among elassifiers - most 
elassifiers agree on the two extremes of the visual interaetion 
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Fig. 2. A montage of randomly chosen inactive galaxies at 1 .5 < z < 2.5 arranged in the fiducial bins of the visual Interaction Metric 
(IM). Each row of five galaxies are randomly selected from all objects in one of the three bins in IM (see Section 2.1.2). A redshift 
is written at the upper left of each panel and IM is written at the lower left. All three color (iJH) images are 6” on a side and are 
identically scaled with a sinh“^ stretch. 


sequence. Intermediate values of IM are more uncertain. For 
example, more than 40% of galaxies with 0.25 < IM < 0.5 
have mean scatters that span across most of the IM sequence 
(0.50 <scatter< 0.75). We are able to trust that a final classifica- 
tion for a galaxy as non-interacting represents a fair consensus 
among visual classifiers, but, if a galaxy shows signs of interac- 
tion, there may be some uncertainty as to the degree. 

For this reason, we bin galaxies by interaction class quite 
coarsely: 0.0 < IM < 0.2, 0.2 < IM < 0.5 and 0.5 < IM < 1.0, 
which represent “isolated”, “interacting” and “merging” objects. 
Figure 2 shows some examples of massive galaxies at 1.5 < z < 
2.5: each row contains randomly chosen galaxies classified into 
the three coarse bins listed above. 

In addition to IM, we identify a refined subset of galaxies 
with interacting or merging activity, which comes from a com- 
bination of visually identified features derived from the visual 
classification catalogs. It consists of galaxies explicitly listed 
as visual mergers or direct interactions (disregarding multiple 
galaxy systems with no clear signs of interactions) by at least 
60% of classifiers, or in which at least 60% of classifiers no- 
ticed tidal tails or double nuclei. This multi-feature subset may 
exclude some very young stage or distant interactions where the 
galaxies have not yet been substantially disturbed, but will be a 
cleaner distillation of actual systems in a broad array of interac- 
tion states, less affected by scatter than the IM. 

We also employ a “dumpiness metric” (CM). CANDELS 
visual classifiers identified clumpy structure by examining im- 
ages of galaxies across multiple bands and were given a choice 
of three alternatives in the GUI. As for IM, we assign numeri- 


cal values to these choices. Objects with smooth structure lack- 
ing any clumps are assigned CM = 0. Those with one or two 
clumps (disregarding any central bulge) are assigned CM = 0.5, 
while objects with three or more clear clumps are given CM = 1 . 
Intermediate values come from averaging the decisions of mul- 
tiple reliable classifiers. 

In the right panel of Figure 1, we show the distributions of 
the scatter in CM calculated in the same manner as as described 
above for IM. The fidelity of the dumpiness measure of galaxies 
is generally better than the interaction measure. Most classifiers 
agree on whether a galaxy is smooth, but there is a small amount 
of classifier scatter on the actual level of dumpiness in galaxies. 
Figure 3 shows some examples of massive galaxies at 1.5 < z < 
2.5: each row contains randomly chosen galaxies classified into 
three coarse bins in CM (0.0 < CM < 0.2, 0.2 < CM < 0.6 and 
0.6 <IM< 1.0). 

In addition to IM and CM, we also use flags set by visual 
classifiers to exclude from our analysis sources with troublesome 
image artifacts and, as needed, those with signs of visual point- 
source contamination. 

It is worth noting that the various measures of structure used 
here are not independent of each other. Some dependencies are 
physical: high Sersic galaxies are generally massive, elliptical 
galaxies which tend to be smooth and, even if interacting, fre- 
quently show only weak signatures of obvious disturbance such 
as tidal tails. Some dependencies are based on the subjectivity 
of visual classification: if a galaxy is judged to be close enough 
to another in the HST images to satisfy a non-zero IM value, a 
disky/clumpy system frequently appears to be more visually dis- 
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Fig. 3. A montage of randomly chosen inactive galaxies at 1.5 < z < 2.5 arranged in the fiducial bins of the visual Clumpiness 
Metric (CM). Each row of five galaxies are randomly selected from all objects in one of the three bins in CM (see Section 2.1.2). A 
redshift is written at the upper left of each panel and IM is written at the lower left. All three color (iJH) images are 6” on a side and 
are identically scaled with a sinh“^ stretch. 


turbed and may be assigned a higher IM than a bulgy/smooth 
system, even if the dumpiness could be related more to the 
gas fraction or inflow history of the galaxy rather than its cur- 
rent interacting state. For example, one of the high IM galax- 
ies in Figure 2 is also found in the high CM row of Figure 3. 
When we compare distributions of individual structural param- 
eters in course of this paper, the reader should keep in mind 
that there may be covariances between the various distributions. 
Disentangling these dependencies is neither trivial nor fruitful, 
since the principal limitation in the determination of trends in 
this work stems from the modest size of the AGN sample. 

2.2. Herschel/PACS Imaging and Photometry 

Our far- infrared data are composed of maps at 70 jum, 100 jum 
and 160 jum from a combination of two large Herschel/PACS 
programs: the PACS Evolutionary Probe (PEP), a guaran- 
teed time program (Lutz et al. 2011) and the GOODS-Herschel 
key program (Elbaz et al. 2011). The combined PEP+GH 
(PEP/GOOD S-Herschel) reductions are described in detail in 
Magnelli et al. (2012). While data at 100 and 160 fim are avail- 
able in both fields, an additional deep map at 70 /urn is also 
available in GOODS-S. The PACS 160, 100 and 70 //m fluxes 
were extracted using sources from archival deep Spitzer MIPS 
24 jum catalogs as priors, following the method described in 
Magnelli et al. (2009); see also Lutz et al. (2011) for more de- 
tails. 3cr depths are 0.90/0.54/1.29 mJy at 70/100/160 jum in 
the central region of GOODS-S and 0.93/2.04 mJy at 100/160 
fim in GOODS-N. The GOODS-S maps are ~ 80% deeper than 


the GOODS-N maps and probe further down the FIR luminosity 
function at all redshifts (Magnelli et al. 2012). 

For practical purposes, we use the monochromatic luminos- 
ity of a galaxy at 60 jum rest (L^o) as a measure of its FIR 
luminosity. The PACS bands cover the rest-frame 60 jum over 
much of the redshift range probed in this work and we estimate 
L^o from a simple log-linear interpolation of PACS measure- 
ments in bands that bracket 60 jum in the rest-frame. The use of 
Leo obviates the need to apply an uncertain correction between 
monochromatic and total FIR luminosities. As a rough guide for 
the reader, a star-forming galaxy with an IR SED similar to M82 
with Leo= 10"^^ L© has a SFR of approximately 80 M©/yr. The 
exact transformation depends on the SED shape, the SF history 
of the galaxy and on many other factors relating to the distribu- 
tion of dust and stars in a system. 

In cases where a mean Leo is desired for a sample con- 
sisting of a mix of PACS detected and undetected sources, we 
follow a technique developed in earlier works from our team 
(Shao et al. 2010; Santini et al. 2012) and briefly outlined here. 
We bin sources in this study by redshift and structure. We em- 
ploy fiducial redshift bins: 0.5 < z < 1.0, 1.0 < z < 1.5, 
1.5 < z < 2.0 and 2.0 < z < 2.5 (though, in some analyses, 
we combine the last two redshift bins to improve statistics). A 
fraction of sources in each bin are detected in two or more PACS 
bands. Leo is calculated for these using their individual redshifts 
and a log-linear interpolation of PACS fluxes. Of the remaining 
sources, some are detected in only one PACS band, while the 
majority are undetected in the FIR data. For the latter, stacks 
were performed at the optical positions of the sources on PACS 
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Table 1. Total AGN / Inactive sample sizes after the application of different selections 


Selection 

0.5 <z< 1.0 

1.0 <z< 1.5 

1.5 <z<2.0 

2.0<z<2.5 

CANDELS/both flelds + M,> 10‘“ Mo (parent) 
parent + good GALFIT 

91/974 

87/919 

89/877 

82/841 

69 / 625 
60 / 570 

34 / 594 
22/460 

CANDELS/GOODS-S + M.> 10'“ Mq (parent) 
parent + visual classification 

52 / 320 
44/293 

46/362 

39/332 

39/211 

33/198 

20/246 

17/229 


residual maps in all available bands, from which mean fluxes 
were measured using PSF photometry. The stacked fluxes in a 
band were averaged with the fluxes of sources only detected in 
that band, weighting by the number of sources in each category. 
This gives mean fluxes for the partially detected and undetected 
AGNs in both bands, from which a mean ^50 was derived using 
the median redshift of these sources. The flnal 60 //m luminosity 
in each bin was computed by averaging over the linear luminosi- 
ties of detections and non-detections, weighted by the number of 
sources. This procedure was only performed for bins with more 
than 3 sources in total. 

Errors on the infrared luminosity are obtained by bootstrap- 
ping. A set of sources equal to the number of sources per 
bin is randomly chosen 100 times among detections and non- 
detections (allowing repetitions), and ^50 is computed per each 
iteration. The standard deviation of the obtained ^60 values gives 
the error on the average 60 jjm luminosity in each bin. The error 
bars thus account for both measurement errors and the scatter in 
the population distribution. 

2.3. Chandra Deep Field (CDF) X-ray catalogs 

Cospatial with the two GOODS survey flelds (Giavalisco et al. 
2004), the Chandra Deep Fields (CDFs) are the deepest pencil- 
beam X-ray surveys in the sky. The 2 Msec exposure in the full 
CDF-North (CDF-N) has produced a point source catalog con- 
sisting of 503 sources (Alexander et al. 2003), while in GOODS- 
South, the recent 4 Msec CDF-South (CDF-S) catalog comprises 
740 sources (Xue et al. 2011). However, only ~ 60% of the full 
CDFs are imaged by CANDELS. 

We have extensively characterized the data and catalogs in 
both flelds, in which careful associations have been made with 
optical and near-IR counterparts, using, where possible, proba- 
bilistic crossmatching models (Luo et al. 2010; Xue et al. 2011). 
In addition to the deep X-ray data, the wealth of deep spec- 
troscopy and multi-wavelength photometric data in the GOODS 
flelds have enabled accurate spectroscopic or AGN-optimised 
photometric redshifts to be determined for the majority of the 
X-ray sources (e.g., Szokoly et al. 2004; Luo et al. 2010). We es- 
timate absorption-corrected hard-band X-ray luminosities (Lx) 
of sources with redshifts using spectral modeling techniques 
(Bauer et al. 2004). As a result of the small area and great 
depth of the CDF exposures, most X-ray sources are low or 
moderate luminosity AGNs - only ~ 5% of the sources have 
logLx(2-10 keV) > 44 erg s“^. These may be thought of as X- 
ray selected equivalents of the local Seyfert galaxy population. 
In this work, we only consider sources with log Lx > 42 erg 
s"\ to prevent contamination from powerful starbursts, in which 
emission from high-mass X-ray binaries can potentially over- 
power the emission from nuclear activity. 

2.4. Multi-wavelength photometry and stellar masses 

We employ multiwavelength galaxy catalogs in both GOODS 
flelds to deflne a general galaxy sample, the properties of which 


we will compare to the AGNs. In GOODS-S, we use the updated 
GOODS-MUSIC database (Santini et al. 2009; Grazian et al. 
2006), while in GOODS-N we use a catalog developed for the 
PEP team using a similar methodology (Berta et al. 2010, 201 1) 
The former catalog selects galaxies with observed magnitudes in 
the HST F850LP band < 26 or in the ISAAC X, band < 23.5, 
while the latter is primarily selected to have K < 24.2. In or- 
der to exclude a surfeit of faint sources with inaccurately red 
colors and masses, we apply an additional cut of F850LP < 26 
in the GOODS-N catalog. For galaxies with no current spectro- 
scopic redshifts, photometric redshifts were determined by fit- 
ting multiwavelength photometry using PEGASE 2.0 templates 
(Fioc & Rocca-Volmerange 1997) in GOODS-S or using the 
EAZY code (Brammer et al. 2008) in GOODS-N. For details on 
the catalog preparation, characterization and photometric red- 
shift estimation, we refer the reader to Santini et al. (2009) for 
GOODS-S and Berta et al. (2010) for GOODS-N. 

We have developed a custom technique to estimate the stel- 
lar masses (ML) in AGNs, by linearly combining galaxy pop- 
ulation synthesis model templates and AGN SED templates to 
fit multiwavelength photometry. For inactive galaxies, we per- 
form a minimization of Bruzual & Chariot (2003) synthetic 
models, assuming a Salpeter Initial Mass Function (IMF) and 
parameterizing the star formation histories as exponentially de- 
clining functions. For AGNs, we also include an AGN template 
from Silva et al. (2004), which accounts for a variable fraction of 
the total light of the galaxy. The AGN template reflects the clas- 
sification of the X-ray source, derived from information about 
its SED and spectrum, where available. For sources classified as 
Type I (broad lines in the spectrum, clear AGN contribution in 
the rest-frame optical and UV), an unobscured AGN SED was 
used. For the rest, a Type II template was used if the estimated 
X-ray absorption column X7/ < 10^"^ cm“^, and a Compton-thick 
template for more heavily absorbed systems, though in practice, 
the choice of the latter two templates makes little difference. For 
further details of the method, performance evaluations, tests and 
limitations, we refer the reader to Santini et al. (2012). 

While AGN are selected by their X-ray emission, we deflne 
our ‘inactive’ galaxy population as all galaxies that are unde- 
tected in X-rays (excluding even those which have log Lx < 42 
erg s“^) and Spitzer/IRAC colors that are unlike those of bright 
AGNs following the criteria of Donley et al. (2012). In practice, 
only a very small fraction of the general galaxy population are 
rejected on the basis of these criteria, but they tend to be mas- 
sive galaxies and could potentially sway the statistics of SF com- 
parisons among such systems by an inordinate degree. We also 
impose a minimum mass of M^= 10^^ M© on both AGNs and 
inactive galaxies. Very few AGNs in our redshift range of in- 
terest lie at lower masses and the GOODS galaxy catalogs be- 
come increasingly incomplete below this mass limit at z ~ 2 
(Santini et al. 2012). 

Galaxies from the multiwavelength catalogs were matched 
to Herschel sources through the positions of the 24 pm priors, 
which are, in turn, tied to IRAC catalogs in the GOODS flelds 
(Magnelli et al. 2009). Therefore, the cross-matching tolerances 
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Fig. 4. Top panels: Physical half-light semi-major axis radii plotted against stellar mass for AGNs (large red points), inactive 
control galaxies (large black points) and the overall galaxy population (small grey points), in four distinct redshift bins. Samples 
from both GOODS fields are combined and one hundred control sample draws are used in all panels of this Figure. The lines in the 
upper panels show parameterisations of galaxy size-mass relationships across redshift from van der Wei et al. (2014); dashed lines 
for late-type galaxies and solid lines for early type galaxies. Vertical histograms compare the size distributions of the AGNs (red 
lines) and the inactive control. Middle panels: Sersic index n plotted against stellar mass M, for AGNs (large red points), inactive 
control galaxies (large black points) and the overall galaxy population (small grey points), in four distinct redshift bins. Samples 
from both GOODS fields are combined here. Vertical histograms compare the size distributions of the AGNs (red lines) and the 
inactive control. In each panel in both rows, the Kolmogorov- Smirnoff (KS) probability {Pks) that the vertical distributions come 
from the same parent sample is shown in the upper left comer. Bottom panels: A comparison of the M, distributions of AGNs 
(red histograms), the inactive control sample (black histograms) and the overall galaxy population (grey filled histograms). The 
importance of stellar mass-matching in constmcting a valid comparison sample is clear, and the performance of our mass-matching 
procedure may be appreciated at a glance. 


are ~ 1”. The rate of spurious crossmatches is ~ 3%, estimated 
from the asymptotic behaviour of the cross-match offsets be- 
tween the I^C catalog and the CANDELS H-band catalog in 
GOODS-S. 


In Table 1, we list the numbers of AGNs and inactive galax- 
ies from the parent samples - sources that lie in the overlap of 
the CANDELS imaging and Herschel/PACS maps, additionally 
restricted by the footprint of the CDFs (for AGNs) or the mul- 
tiwavelength catalogs (for the inactive galaxies). Along with the 
parent sample numbers, we also list numbers of AGNs and in- 
active galaxies that satisfy our criteria for good GALFIT mea- 
surements (in both GOODS fields) and good visual structure es- 
timates (in GOODS-S only). 


2.4.1. Bootstrapped stellar mass-matched control samples 

Moderately luminous AGNs are inherently rare among galaxies. 
In the small fields considered in our study, only few to several 
tens of AGNs that satisfy our the minimum criteria for a valid 
structural measurement may be found in each redshift bin. Direct 
constraints on the relationships between SFR, redshift and struc- 
ture based solely on the AGNs themselves are severely limited 
by small number statistics. 

On the other hand, the inactive galaxy population far out- 
number X-ray selected AGNs in both fields and at all redshifts. 
We take advantage of this by creating multiple control samples 
of inactive galaxies matched to the AGNs using Monte-Carlo 
bootstrap techniques, from which we constrain the measure- 
ments and distributions of structural and SF properties more rig- 
orously, as well as account for the uncertainties arising from the 
small sample size of the AGNs and complex scatter associated 
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Fig. 5. Distributions of the light profile parameters from GALFIT fits for AGNs and inaetive eontrol galaxies. The three rows of 
panels show histograms in Sersie index (n), ellipse axis ratio (q) and physieal semi-major axis radius (Rg) from top to bottom 
respeetively. Panels left to right span four distinet bins in redshift. Samples from both GOODS fields are eombined. 1000 draws 
of a mass-matehed eontrol sample are analysed to determine the lcr/2cr seatter in the distributions of the parameters for inaetive 
galaxies, shown as dark/light grey zones in the histograms. Red open histograms show the parameter distributions for the AGNs. 


with the struetural measurements. The observed measurements 
and distributions shown by the AGNs may then be eompared to 
those of the eontrol samples. Sinee the uneertainties in the latter 
are also estimated, we ean aseertain whether or not the observed 
AGNs are eonsistent with being drawn from the inaetive galaxy 
population. 

We mateh galaxies to AGNs on the basis of stellar mass. 
All AGNs and inaetive galaxies in eaeh of the fidueial redshift 
bins are further binned into narrow mass intervals, of AM=^= 0.2 
dex for M^< 10^^-^ M© and AM^= 0.5 dex for M^> 10^^-^ M©. 
The inerease of the matehing toleranee at high M, is due to the 
paueity of high mass galaxies, as well as the high AGN inei- 
denee at these masses, whieh reduees the pool of eontrol galax- 
ies substantially. For eaeh AGN in a redshift and mass bin, we 
randomly ehoose one eounterpart from the eorresponding set of 


inaetive galaxies, allowing for duplieations. This yields a single 
set of inaetive eontrol galaxies of the same number as the AGNs, 
sharing their stellar mass distribution within the matehing tol- 
eranee. We repeat this proeess hundreds of times to get multi- 
ple independent eontrol samples. These are used to determine a 
“statistieal distribution” of any parameter of interest, whieh en- 
eapsulates both the mean distribution and the statistieal seatter in 
the distribution eoming from real and sampling varianee of both 
the inaetive galaxy population and the measurements of the pa- 
rameter itself. We then evaluate whether the distribution of the 
parameter for the AGNs is eonsistent with arising as a single 
draw of the eontrol sample. In this sense, the approaeh taken in 
this work is purely eomparative and no attempt is made to eorreet 
the distributions for ineompleteness or biases assoeiated with the 
photometrie euts and the quality euts applied for struetural mea- 
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surements. Since both AGNs and galaxies are tied to the same 
multi-wavelength catalogs, the same biases are expected to ap- 
ply to both populations. 

Throughout the rest of this paper, we will use the term “in- 
active galaxies” to specifically denote the population of massive 
galaxies without detectable nuclear activity with a stellar mass 
distribution matching those of the AGNs, i.e., the stellar mass- 
matched control samples described above. 

3. Setting the stage: Structural patterns for AGNs 
and inactive galaxies 

We start by highlighting a few important structural relationships 
inherent for galaxies, as uncovered by our morphological mea- 
surements. These serve as context for understanding differences 
between AGNs and inactive galaxies in terms of their relation- 
ships between SF and galaxy structure. 

3. 1. Size-mass relationship 

Galaxies show a clear trend towards larger sizes at higher stel- 
lar mass, a relationship that is discernible to z > 2. We consider 
the size-mass relationship in the upper panels of Figure 4, where 
we plot H-band half-light radii Rg vs. M, for galaxies and AGNs 
combining the two GOODS fields. As a guide to the eye, lines 
are used to illustrate the typical size-mass relationships for star- 
forming galaxies and quiescent galaxies as recently determined 
from the CANDELS-based study of vander Wei et al. (2014). 
Since current star-formation is correlated with the light/mass 
profile of galaxies, early- and late-type galaxy populations with 
different Sersic indices (n) are also differentiated by their char- 
acteristic size-mass relationships. High n galaxies show a much 
steeper relationship between size and mass and are generally 
smaller than low n galaxies at the same M^. In addition, the size- 
mass relationship of high n galaxies evolves more strongly to 
z - 2 (Trujillo et al. 2006, 2007; Bruce etal. 2012). The dif- 
ferences between these relationships and their evolution place 
important constraints on the buildup of stellar mass in galaxies 
fromz = 2. 

The AGNs and 100 sets of mass-matched control galaxies 
are shown in this diagram as larger red and black points respec- 
tively. For the most part, AGNs occupy the same range and scat- 
ter in size as shown by inactive galaxies of the same stellar mass. 
This is shown more clearly through the vertical histograms on 
the side of each panel, which compare the size distributions of 
the AGNs and the control galaxies. The similarity of the distribu- 
tions is measured using a two-sided Kolmogorov- Smirnoff (KS) 
test and the probability that the two distributions are drawn from 
the same parent population (Pks) is listed in the upper left comer 
of each panel. We see some systematic variation in the relative 
size distributions with redshift. The KS tests suggest that the two 
distributions are very consistent between z = 1 toz = 2.5, but 
differs in the lowest. At 0.5 < z < 1.0, AGN hosts are larger 
than the control sample. However, a detailed look in Section 3.3 
suggests only minor differences, which may be driven by cosmic 
variance. 

3.2. Sersic index-mass relationship 

The presence of two size-mass relationships among galaxies is 
related to differences in their light profiles, as parameterised by 
the Sersic index (n) (Trujillo et al. 2006) This is believed to be 
driven by their different dynamical evolution (e.g. Baugh et al. 


1996; Mo et al. 1998; Naab et al. 2009). Galaxies with n ^ 1 
have primarily exponential light profiles and are likely to be 
disk-dominated. Since disks are dynamically colder systems 
with substantial rotation, low Sersic galaxies are unlikely to have 
suffered a violent event within several dynamical times of our 
view of them. On the other hand, galaxies with ^ ^ 4 are clas- 
sical spheroids, which are dynamically hot and have low rota- 
tion. This suggests that their most critical recent evolutionary 
event was violent enough to redistribute any cold components 
that may have existed in the progenitors of these galaxies. This, 
and the relative importance of minor mergers and recent gas ac- 
cretion, play a role in the interpretation of the different size-mass 
relationships shown by low and high Sersic galaxies. Some com- 
plexity to this picture comes from suggestions that modest Sersic 
indices could be achieved by disk-dominated galaxies in the high 
redshift Universe (e.g. van der Wei et al. 2011). 

In the middle panels of Figure 4, we plot n vs. M, for 
galaxies and AGNs for the two GOODS fields combined. At 
M^< 10^^ M© (with possibly some differences over redshift), 
the vast majority of galaxies have low Sersic indices ofn< 2.5. 
Towards M^> 10^^ M©, the typical Sersic index of galaxies in- 
creases considerably, reflecting the increased spheroidal domi- 
nance of high mass galaxies seen both in the local and distant 
Universe (e.g., Buitrago et al. 2013). There is clear change in 
the fraction of high n galaxies with redshift over the range we 
study here. At 0.5 < z < 1.0 and M^> 10^^ M©, 35% of galax- 
ies have n > 2.5, which changes to 18% among equally massive 
galaxies at2.0<z<2.5. This is a well-established result in the 
CANDELS era (e.g. van der Wei et al. 201 1). 

As before, we show the AGNs and control galaxies using 
larger red and black points, and compare their Sersic index dis- 
tributions visually using vertical histograms in each panel and 
analytically with two-sided KS tests. At z < 1.5, AGNs and in- 
active galaxies have statistically indistinguishable distributions. 
However, at higher redshifts, AGNs show a significant surfeit of 
high Sersic systems. In the next subsection, we examine the dif- 
ferences between AGNs and inactive galaxies more carefully to 
gain some insight into the causes for these observed differences. 

3.3. Light profile parameter distributions 

The large number of inactive galaxies from the parent catalogs 
enables the construction of statistical distributions of n, Rg and 
the ellipse axis ratio q for the mass-matched inactive control 
sample, combining both GOODS fields to maximise the sam- 
ple size and minimize cosmic variance. In Figure 5, each row 
shows histogram distributions for a different light profile param- 
eter output by GALFIT. In each panel, the dark/light grey regions 
show, for each histogram bin, the Icr/lcr range in the number of 
inactive galaxies in that bin, determined from 1000 control sam- 
ple draws. The distribution for the AGNs are shown as red open 
histograms. If the AGN histogram lies within the dark grey re- 
gions, their distributions are completely consistent with those of 
inactive galaxies within 1 standard deviation. 

By and large, the distributions of all the structural parameters 
among AGNs are not radically dissimilar from the control sam- 
ple across all redshifts. Closer scrutiny suggests some systematic 
differences. The distributions of Sersic index of the AGNs, are 
broadly consistent with those of the control galaxies at z < 1.5 
but deviate significantly towards higher n at z > 1.5. At these 
high redshifts, the shapes of the AGNs are a bit rounder, with 
a q distribution skewed slightly towards higher values than the 
inactive galaxies. AGNs also tend to be a bit larger than inactive 
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Fig. 6. GALFIT-derived Sersic index vs. absorption-corrected 2-10 keV X-ray luminosity of AGNs from both Chandra Deep Fields 
combined (white circle points). The contours show the density of points from 1000 draws of a mass-matched control sample of 
inactive galaxies. To place the control galaxies on this plot, the X-ray luminosity assigned to each control galaxy is that of the AGN 
to which it is matched, producing, in effect, a large “simulated AGN” sample with the same mass and X-ray luminosity distributions 
as the real AGNs. The four contour levels encompass 90%, 75%, 50% and 25% of the simulated AGNs in this plane. 
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Fig. 8. A comparison of a typical galaxy spectrum (red solid 
line) and AGN spectra with different levels of extinction (blue 
lines). The base galaxy spectrum is the ‘Sa’ template from the 
SWIRE library (Polletta et al. 2007) scaled to correspond to 
H = 24 magnitude if placed at z = 2. The base AGN spec- 
trum (blue solid line) is stitched together from composite radio- 
quiet QSO spectra in the UV from Telfer et al. (2002) and the 
optical from Vanden Berk et al. (2001), scaled to correspond 
to Lx= 10"^^-^ erg s“^ using optical/X-ray relation defined by 
Lusso et al. (2010); the error bars show the expected scatter in 
the relation. The blue dashed lines show the effects of Calzetti- 
law extinction on the base AGN spectrum with E(B-V) as indi- 
cated on the right below the corresponding extinguished spec- 
trum. The approximate position of the observed F850LP (z) and 
F160W (H) bands at z = 2, relevant for the interpretation of 
Figure 6, are shown as grey bands. 


galaxies at 0.5 < z < 1.0 and a bit smaller than inactive galaxies 
at 1.5 < z < 2.5. 

An interesting trend to note is that the light profile distri- 
butions of AGNs change very little across redshift, unlike that 
of the inactive galaxies. Between z = 0.5 and z = 2.5, the 


Sersic index distribution of the AGNs has a roughly fixed me- 
dian value of ~ 2.7, while that of the inactive galaxies ranges 
from 3.0 at 0.5 < z < 1.0 to ~ 1.5 at 2.0 < z < 2.5. In 
parallel, inactive galaxies also become more elongated, a sign of 
the increased prominence of massive disks with redshift. AGNs, 
however, remain rounder, with only a slight shift towards more 
elongated profiles between the low and high redshift bin. The 
median sizes of the control galaxies increases slightly with red- 
shift, while those of the AGNs decreases slightly with redshift 

Taken together, these systematics can be attributed to two 
possible effects. A small amount of point source contamination 
from the nucleus may be affecting AGN light profiles, mostly in 
the high redshift bin where the AGNs are more luminous and 
the H-band traces the rest-frame B and V bands (4400-6000 A). 
Alternatively, there may be a preference for AGNs to be found 
in more bulge-dominated hosts, as suggested by earlier visual 
classification studies (Kocevski et al. 2012). Both effects lead to 
a greater concentration of light coming from the central region 
of the galaxy, resulting in the small differences in light profile 
structure we observe. 

3.3.1. Possible origins of the central light excess 

An AGN can be a profuse source of energy emitting across most 
wavelengths. Unextinguished UV-to-optical continuum emis- 
sion from a nucleus originates from the hot accretion disk and is 
very blue, rising in fiux rapidly at shorter wavelengths towards 
a putative peak in the extreme UV. On this continuum may be 
superimposed emission lines from the various AGN-ionized re- 
gions. Regarding its effect on the appearance of galaxy structure, 
unextinguished AGN appear as nuclear point sources with blue 
colors. However, dust along the line of sight to the nucleus can 
redden the point source, making it harder to distinguish against 
the light from the bulge of the host galaxy. 

The infiuence of strong AGN point sources on galaxy struc- 
tural measurements can be quite profound. Even point source 
fractions as low as a few tens of percent can lead to appre- 
ciably bulgier, rounder and more centrally concentrated effec- 
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Fig. 7. Distributions of the Sersic index {n) from GALFIT fits for AGNs and inaetive eontrol galaxies in two HST bands: the 
ACS/F850LP band (bottom panels) and the WFC3/F160W band (top panels). Panels left to right span three distinet bins in redshift. 
Only the subset of sourees eontained in the FIREWORKS eatalog are here. 1000 draws of a mass-matehed eontrol sample are 
analysed to determine the lcr/2cr seatter in the n distributions for inaetive galaxies, shown as dark/light grey zones in the histograms. 
Red open histograms show the parameter distributions for the AGNs. Sersie indiees are typieally lower in the bluer F850LP band. 
Espeeially at 1.5 < z < 2.5, the n histograms of the AGNs are more eonsistent with those of the inaetive galaxies from the F850LP 
fits. 


tive light profiles (Pieree et al. 2010). It is worth noting, given 
the nuelear luminosities of the AGNs in the GOODS fields 
(as estimated from their X-ray emission), only a small fraetion 
should be powerful enough to severely eontaminate their host 
galaxy light, even if unextinguished by dust (see the Appendix 
of Rosario et al. (2013 a) for a diseussion). Earlier GALFIT stud- 
ies have revealed frequent red nuelear exeesses among AGNs in 
the GOODS field (Simmons et al. 2011; Sehawinski et al. 201 1), 
attributable to obseured AGN point sourees. These sourees are 
generally mueh fainter than their hosts and would not strongly 
alter their measured strueture. However, pervasive weak nuelear 
exeesses among AGNs eould explain some of the mild system- 
atie trends we find in the light profile parameter distributions of 
AGN hosts. 

In Appendix A, we undertake a brief study of the evidenee 
for nuelear eontamination of the light profiles of AGNs at 1.5 < 
z < 2.5. We fit AGNs and mass-matehed inaetive galaxies us- 
ing a eombination of a Sersie galaxy model and a eentral PSF. 
The inelusion of a point souree in the GALFIT fits leads to 
a lower effeetive Sersie index distribution for the AGN hosts 
(Figure A.l), but it also systematieally lowers the effeetive in- 
diees for inaetive galaxies. As strong nuelear emission is not ex- 
peeted in the latter population, we eonelude that sueh two eom- 


ponent fits are subtraeting away the light of a eentral bulge as 
well. This is not unreasonable given the similar sizes of bulges 
and WFC3/F160W PSFs at z ~ 2 (e.g. Bruee et al. 2012). An 
aeeurate treatment of the influenee of AGN nuelear emission re- 
quires multi-eomponent fits, as well as eareful simulations with 
the addition of fake point sourees to real images of massive in- 
aetive galaxies in order to understand the systematies inherent in 
the eonsideration of a seeond struetural eomponent in GALFIT. 
This will be pursued in future work built upon this study. We 
ean, nonetheless, perform some simple tests for the importanee 
of strong blue exeesses to the light profiles of AGNs. 

In Figure 6, we plot the Sersie indiees of AGNs against their 
absorption-eorreeted 2-10 keV X-ray luminosity Lx, eombining 
the samples from both GOODS fields. The underlying eontours 
show the distribution in this plane of 1000 draws of a mass- 
matehed eontrol sample of inaetive galaxies. To be plaeed in this 
Figure, eaeh eontrol galaxy is assigned the X-ray luminosity of 
the AGN to whieh it is matehed, generating, in effeet, a large 
“simulated AGN” sample. If the light profiles of AGNs were 
independent of nuelear luminosity and identieal to those of in- 
aetive galaxies, they are expeeted to seatter in this plane in the 
same fashion as the eontrol sample. This is the ease at z < 1. At 
higher redshifts, there is a definite tendeney for the more lumi- 
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nous AGNs to occur in systems with higher n. In particular, one 
may note that AGN hosts with n >4, which may be seen at z < 1 
across Lx, are more common in sources with Lx> 10^^ erg s"^ 
in the two higher redshift samples. This Figure seems to suggest 
that the more luminous AGN hosts at z > 1 .5 are associated with 
high Sersic systems, perhaps due to higher levels of AGN con- 
tamination, or because such luminous AGNs are preferentially 
found in hosts with larger spheroids. A similar result was found 
in the visual study of X-ray AGN hosts in CANDELS/GOODS- 
S by Kocevski et al. (2012). 

We perform another simple test of these alternatives by em- 
ploying GALFIT fits of the same galaxies in a bluer band. If blue 
nuclear contamination is the reason for the systematic deviation 
of n in AGN hosts at z > 2, and assuming that AGN spectra 
do not change systematically with redshift, then we should see a 
qualitatively similar deviation at 0.5 <z< 1.0 when we examine 
n distributions from F850LP (z-band) fits, since this band traces 
the re St- frame B-band for these galaxies. In addition, the devi- 
ation should be much more pronounced for the higher redshift 
systems where F850LP traces rest-frame UV light. 

In Figure 7, we plot the n distributions of galaxies and AGNs 
in the GOODS-S/FIREWORKS catalog, for which consistent 
fits using the same setup are available in both H- and z-bands 
(see Section 2.1.1). A quick inspection of the top row of panels 
in this Figure and in Figure 5 will show that the basic H-band n 
distributions and trends are qualitatively preserved whether one 
is considering the complete AGN sample or just this subset in 
GOODS-S. Comparing the two rows of panels in Figure 7, we 
see that, while the overall distributions for both AGNs and con- 
trol are slightly diskier in the z-band, the AGNs do not show 
more concentrated light profiles in these bands, as one would 
have expected from prominent blue central point sources. In 
particular, the AGNs appear to be diskier than the control at 
0.5 < z < 1.0, not bulgier. Also both AGNs and inactive galaxies 
show strongly disky light profiles at z > 1.5, arguably closer to 
that of the control in the rest- frame UV (z-band) than in the rest- 
frame optical (H-band). The greater similarity of the AGN and 
control z-band n distributions at z > 1.5 suggests that the excess 
light in the AGN hosts is at least as red or redder than the typical 
colour of the control galaxies, since it only becomes prominent 
at optical wavelengths. Therefore, if nuclear light is the major 
source of the central excess, it must be reddened considerably 
by dust along the line-of-sight to the nucleus. However, this will 
also act to preferentially extinguish the light from the nucleus, 
decreasing the contrast of nuclear emission with respect to the 
host galaxy. 

We explore the consistency of the notion that wide-spread 
reddened AGN light can infiuence the centres of the more lu- 
minous host galaxies at z ~ 2 using Figure 8, in which we 
contrast a typical galaxy spectrum to AGN spectra with vary- 
ing degrees of extinction. The galaxy template shown in the 
Figure is a typical ‘Sa’ galaxy spectral type from the SWIRE 
template library (Polletta et al. 2007), normalised to a luminos- 
ity at z = 2 that corresponds to an observed frame magnitude of 
H = 24. The AGN spectrum comes from a composite HST/UV 
spectrum of radio-quiet QSOs (Telfer et al. 2002) stitched to a 
composite optical spectrum from the Sloan Digital Sky Survey 
(Vanden Berk et al. 2001) (the tabulated spectrum is available at 
http://www.pha.jhu.edu/~rtl9/composite/). Since it is empirical, 
this base spectrum already includes a mild degree of extinction. 
The QSO spectrum is normalised to the characteristic luminosity 
of an AGN with Lx= 10"^^-^ erg s“^, using the X-ray to optical 
relation for Type 1 AGNs from (Eusso et al. 2010). Different lev- 
els of extinction applied to the AGN spectrum yield the spectra 
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Fig. 10. Comparison of the observed fractions of AGNs (red 
points) and inactive galaxies (black points) that are classified 
as interacting/merging using the multi-feature approach, in the 
three fiducial redshift bins (edges shown with dotted lines). The 
inactive galaxy fraction is the median value from 1000 draws 
of a mass-matched control sample to the AGNs in that redshift 
bin, while the error bars show the Icr scatter from those draws. 
Therefore, the AGNs have an enhancement in the observed in- 
teraction fraction over inactive galaxies at z < 1 (by ~ 2cr), but 
are comparable to inactive galaxies at z ~ 2. 

shown with dashed lines. An extinction law from (Calzetti et al. 
1994) was assumed, but the conclusions are not strongly depen- 
dent on this choice, or on the choice of the galaxy template used 
in this exercise. 

An X-ray AGN at z = 2 with the typical luminosity of those 
in our sample can produce roughly as much light at 5000 A 
as a galaxy at the faint limit considered in this structural study 
(H = 24). However, even modest amounts of extinction quickly 
weaken the relative blue light of such an AGN, leaving the light 
from the host galaxy dominant. For an AGN to simultaneously 
account for more than 10% of the total light of a system yet 
be as red as a galaxy spectrum, it would have to be consider- 
ably more luminous than Lx= 10^^ erg s"\ and only a hand- 
ful of such sources are found in our AGN sample. This analysis 
does not preclude that these rather low luminosity systems could 
have a much lower X-ray/optical luminosity ratio (a ox) than cal- 
ibrations based on bright QSOs (Vignali et al. 2003; Eusso et al. 
2010) or a much redder intrinsic spectrum than those of QSOs, 
such as those used to construct the template shown in Figure 
8. Most studies, however, find a higher aox in low-luminosity 
AGN (e.g., Steffen et al. 2006). 

All together, these tests lend some support to the notion that 
high redshift AGNs have a more prominent bulge than equally 
massive inactive galaxies, and that the more luminous AGNs are 
among the bulgiest. We develop the implications of these results 
in Section 5. 

3.4. Distributions of the Interaction Metric 

A comparison of the distributions of the visual interaction met- 
ric (IM) allows us to quantify the fraction of AGNs in interact- 
ing or merging systems relative to inactive galaxies and explore 
the relationship between galaxy interactions and the triggering 
of nuclear activity. In Figure 9, we plot the statistical distribu- 
tions of IM for inactive galaxies, divided coarsely into the three 
bins delineating “isolated”, “interacting” and “merging” systems 
(Section 2.1.2). The histograms come from 1000 draws of mass- 
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Fig. 9. Distributions of the visual Interaction Metric (IM) for AGNs and inactive control galaxies. Panels left to right span three 
distinct bins in redshift. 1000 draws of a mass-matched control sample are analysed to determine the lcr/2cr scatter in the IM 
distributions for inactive galaxies, shown as dark/light grey zones in the histograms. Red open histograms show the parameter 
distributions for the AGNs. Interacting and merging systems are a minor fraction of galaxies and AGNs at all redshifts. AGNs 
are generally equally likely to be in an interacting host as inactive galaxies, though a small excess of AGNs in mergers may be 
discernible at low redshifts. 
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Fig. 11. Distributions of the visual Clumpiness Metric (CM) for AGNs and inactive control galaxies. Panels left to right span three 
distinct bins in redshift. 1000 draws of a mass-matched control sample are analysed to determine the lcr/2cr scatter in the CM 
distributions for inactive galaxies, shown as dark/light grey zones in the histograms. Red open histograms show the parameter 
distributions for the AGNs. Visually clumpy systems are a minor fraction of galaxies and AGNs at all redshifts. At 0.5 < z < 1.0, 
AGNs are significantly more likely to be in clumpy galaxies. By 1.5 < z < 2.5, the CM distributions are similar. 


matched control samples in GOODS-S. These may be compared 
to the IM distributions of the CDF-S AGNs (red histograms) in 
the three fiducial redshift bins. 

A key point to notice is that the majority of sources through- 
out the redshift range 0.5 - 2.5 lie in “isolated” systems, with 
IM= [0.0, 0.2]. The fraction of galaxies classified as “interact- 
ing” or “merging” (IM= [0.2, 1.0]) appears to evolve with red- 
shift (from 18% at 0.5 < z < 1.0 to 30% at 1.5 < z < 2.5, 
estimated directly from the control samples). However, as noted 
before, morphological K-corrections makes such evolution po- 
tentially hard to interpret. 

From Figure 9, we find a significant difference in the inter- 
action properties of AGNs and inactive galaxies only at z < 1. 
In the low redshift bin, there is an excess of “merging” systems 
among AGNs (several cr) at the expense of a smaller number of 
isolated systems. In both the higher redshift bins, the IM distri- 


butions of AGNs and inactive galaxies are consistent within the 
scatter. 

Another approach to test the relevance of mergers is to con- 
sider the fraction of AGNs and inactive galaxies in the ’’multi- 
feature” subset of merging/interacting galaxies, described in 
Section 2.1.2. In Figure 10, we plot as black points the me- 
dian fractions of the mass-matched inactive control sample from 
1000 draws in each of the three redshift bins. The errors on these 
fractions represent the scatter among inactive galaxies from the 
draws. The observed fraction of interactions among the AGNs 
are shown as red points, and may be compared to the median 
and scatter of the interaction fractions among inactive galaxies. 
In parallel to what we find in the IM distributions, the AGNs 
show a significant excess of interacting systems at z < 1 , but 
are consistent within the scatter with the interaction fractions of 
inactive galaxies atz ~ 2. 
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Fig. 12. An illustration useful for understanding Figures 13, 14, 
15. The solid eurve shows a qualitative sehematie SFR distri- 
bution of a population of galaxies, for example, massive galax- 
ies at z ~ 2 where SF galaxies outnumber the quieseent pop- 
ulation. It is assumed that the ^60 distribution is identieal and 
traeks the SFR. The measurable ^60 (or SFR) distribution of 
galaxies individually deteeted in the Hersehel/PACS maps is 
shown as a representative grey histogram. At these redshifts, 
the PEP+GOODS-Hersehel PACS eatalogs are sensitive to FIR 
luminosities that traee just deeper than the ridgeline of the SF 
Sequenee (e.g. Rosario et al. 2013b), but the eompleteness of the 
eatalogs at this limit is low (~ 30%) leading to a shallow eutoff in 
the histogram at the faint end. The mean ^60 of the entire popula- 
tion eomes from staeks into the PACS maps (see Seetion 2.2) and 
is shown as a dashed vertieal line. Sinee staeking is a purely lin- 
ear proeess, the mean is shifted towards the star-forming end in 
this diagram. Additionally, the mean L6o of both PACS deteeted 
and undeteeted galaxies ean be determined, shown as a star and 
eirele point as in the later Figures. The sehematie also reveals 
how the fraetion of galaxies that are FIR-deteeted is elosely re- 
lated to the fraetion of moderate and strongly SF galaxies in the 
population. 

3.5. Distributions of Clumpiness Metric 

A measure of the level of disturbanee in a galaxy is the elumpi- 
ness metrie (CM), whieh is based on the visual prominenee 
of elumps in galaxies. Galaxy simulations suggest that the 
torques driven by large elumps ean inerease the inflow of gas 
onto SMBHs, triggering AGN aetivity (Boumaud et al. 2011; 
Gabor & Bournaud 2013). X-ray staeking studies of elumpy 
galaxies at z ~ 1 support this theoretieal insight (Bournaud et al. 
2012 ). 

From 1000 draws of mass-matehed eontrol samples, we de- 
termine the statistieal distributions of CM for inaetive galaxies 
in our fldueial redshift bins, as shown in Figure 11. There is a 
sharp dropoff in the distribution with CM; most galaxies only 
show low levels of visual elumpiness. The fraetion of elumpy 
systems (defined here as CM > 0.2) inereases mildly with red- 
shift from 25% at 0.5 < z < 1.0 to 32% 1.5 < z < 2.5. This 
is expeeted given the eurrent understanding of galaxy evolu- 
tion where more turbulent high redshift galaxies eontain more 
elumpy disks. However, sinee the bluest band (F606W) used 
for the visual assessment of elumpiness ehanges from rest- frame 
~ 3 500 A in the low redshift bin to ~ 2000 A in the high redshift 
bin, a morphologieal k-eorreetion eould also play some role in 
this apparent evolution. 

The distribution of elumpiness in AGNs is overplotted in the 
Figure as a red solid line. In many ways, the trends found in IM 


in Figure 9 are also refleeted in CM. At 0.5 < z < 1.0, AGNs are 
more likely to be found in elumpy galaxies, at a signifleanee of 
several times the seatter. At these redshifts, we estimate a elumpy 
fraetion among AGNs of 39'^^%. By 1.5 < z < 2.5, the fraetion 
of elumpy AGNs drops to 28;^^%, though it is formally eonsis- 
tent with the elumpy fraetion of AGNs all redshifts. This fraetion 
is eomparable to the elumpy fraetion of inaetive galaxies at the 
same redshifts, and an examination of the right panel in Figure 
1 1 also shows that AGNs and inaetive galaxies have similar CM 
distributions at higher z. This may be eontrasted with their rather 
different distributions at low redshift. Broadly, there appears to 
be a preferenee for the elumpy fraetion of AGNs to ehange less 
slowly with redshift, akin to their light profile parameter distri- 
butions. This implies that the eharaeteristie strueture of AGN 
hosts may remain approximately invariant with redshift between 
z = 2.5 andz = 0.5. 


4. SFR-structure relationships 

Valuable insight into the seenarios that link AGN fueling and 
host galaxy properties eomes from eombining information about 
SFR and host strueture. In this seetion, we explore these re- 
lationships using the deep Hersehel/PACS FIR data in the two 
GOODS fields. We employ two main types of measurements in 
this analysis. Firstly, we use mean FIR luminosities of AGNs 
and inaetive eontrol galaxies. These are estimated direetly from 
PACS fluxes for the subset of sourees deteeted in the PACS 
maps, and from staeks on the residual maps for the PACS- 
undeteeted sourees (Seetion 2.2). Seeondly, we employ the FIR 
deteetion fraetion, defined as the fraetion of sourees deteeted 
in both the 100 and 160 pm PACS maps. The PEP+GOODS- 
Hersehel photometry deteets massive galaxies that lie on or 
just below the SF sequenee out to z = 2.5 (Rosario et al. 
2013b; Magnelli et al. 2013). Therefore, the FIR deteetion frae- 
tion among a eertain population of galaxies is a measure of 
the fraetion of moderate and strongly star-forming galaxies in 
that population. In absolute terms, this measure is sensitive 
to the depth and eompleteness of the PACS eatalogs, and the 
redshift-dependent FIR luminosity funetion. We only use it here 
as eomparative tool to understand differenees between AGNs 
and inaetive galaxies, and speeifleally warn against the over- 
interpretation of any trends with redshift. 

4. 1. SFR as a function of Sersic index 

In Figure 13, we eompare, in the upper row of panels, the mean 
L^o of AGNs (red points) and mass-matehed eontrol galaxies 
(blaek points) in bins of redshift and Sersie index n. The er- 
rors on the mean L^q for AGNs eome from a bootstrap resam- 
pling of the AGN sample itself (Seetion 2.2), while the errors on 
the mean L^q for the inaetive galaxies are determined from 100 
draws of mass-matehed eontrols. 

For the eontrol sample, we And a elear and eonsistent drop in 
the mean Leo with n in the low and intermediate redshift bins, al- 
though the trend appears to flatten out at 1.5 < z <2.5. This is as 
expeeted - most of the SF in the loeal Universe and at higher red- 
shifts is assoeiated with galaxies having substantial disks (e.g. 
Kauffmann et al. 2003; Wuyts et al. 2011; Lee et al. 2013), al- 
though at z ~ 2, a sizeable population of blue star-forming eom- 
paet spheroids is also found (Barro et al. 2013). In eontrast to the 
inaetive galaxies, the AGNs show a mueh flatter dependenee of 
mean Leo with n, sueh that low Sersie galaxies and AGNs have 
eomparable FIR luminosities, but, espeeially at z < 1.5, high 
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Sersic AGNs are significantly more luminous than high Sersic 
galaxies. In essence, this means that the enhancement in mean 
SFR among AGNs is more pronounced in bulgy or early-type 
galaxies. 

In the lower panels, we compare the ^50 distributions of 
AGNs and inactive galaxies, splitting them by redshift and 
coarse bins in Sersic index: 0.5 < w < 2.5 cover disk-like ex- 
ponential light profiles, while 2.5 < n < 5.5 include bulge-like 
De Vaucouleurs profiles. The histograms show the Z 60 distribu- 
tions of PACS-detected sources. As before, the dark/light grey 
shaded regions show the l(r/2cr scatter of the mass-matched con- 
trol sample, and the AGNs are represented by colored open his- 
tograms. The mean ^60 of the PACS-detected objects are shown 
as star points, where error bars represent one standard deviation 
in the mean ^60 of the control sample. The ^60 distributions of 
PACS-detected AGNs are statistically very similar to those of the 
control sample, except among isolated galaxies at high redshift, 
where the AGNs show a slightly lower fraction of FIR bright 
systems. 

We also compare the mean of PACS-undetected sources 
based on stacks, which are shown in Figure 13 as large circle 
points. Again, where measurable, the AGNs have comparable 
mean FIR luminosities to the inactive galaxies, except in the low 
redshift bin, where they appear to be more luminous by a fac- 
tor of a few. Finally, we show the PACS detection fraction of 
AGNs and inactive control galaxies as cross points in the insets 
to each panel of the lower Figure. As has been noted by us in 
earlier studies, AGNs as a whole are significantly more likely to 
be detected in PACS than equally massive inactive galaxies, sug- 
gesting that nuclear activity prefers SF hosts (Santini et al. 2012; 
Rosario et al. 2013b). The biggest consistent difference between 
AGNs and control galaxies can be found in the detection frac- 
tions. For low n, AGNs have a mildly elevated FIR detection 
fraction over the control, comparable to the full AGN sample. 
However, high n AGNs are several times more likely to be de- 
tected in PACS than inactive galaxies with similar structure. This 
is the likely reason for the enhanced mean seen among bulgy 
AGNs in the upper panels. 

4.2. SFR as a function of Interaction Metric 

Our analysis of Section 3.4 suggests that merger rates among 
AGNs at z < 1 are enhanced. If substantial fueling of AGNs 
and bursts of SF are both instigated by gas-rich mergers, nuclear 
activity is expected to be preferentially found in those interact- 
ing systems that also exhibit strong SF. We test this notion in 
GOODS-S using the visual interaction metric to separate iso- 
lated, interacting and merging systems. 

In Figure 14, we compare, in the upper row of panels, the 
mean L^o of AGNs (red points) and mass-matched control galax- 
ies (black points) in bins of redshift and IM. In general, we find 
that the FIR luminosities are not strongly dependent on IM for 
both classes of sources. There is a rise between “isolated” and 
“merging” systems, by a factor of ~ 0.5 dex in the two higher 
redshift bins. The mean L^o of the AGNs is somewhat enhanced 
with respect to inactive galaxies at z < 1.5, by ~ 0.3 dex. This is 
broadly consistent with the level of FIR enhancement found for 
all X-ray AGNs in Santini et al. (2012). There is a hint that the 
enhancement may be higher among “interacting” galaxies in the 
lowest redshift bin, but this pattern is not observed at intermedi- 
ate redshift. At 1.5 < z < 2.5, AGNs and inactive galaxies show 
similar mean Leo. 

In the lower panels, we compare the Leo distributions of 
AGNs and inactive galaxies divided by redshift and into bins of 


0.0 <IM< 0.2 (isolated galaxies) and 0.2 <IM< 1.0 (all interact- 
ing and merging systems for sufficient statistics). The Leo distri- 
butions of PACS-detected AGNs are generally indistinguishable 
from those of the control sample, except among isolated galax- 
ies at high redshift, where the AGNs show a slightly lower frac- 
tion of FIR bright systems. The mean Leo values (star points) 
are consistent with the visual appearance of the distributions and 
also indicate the significance of the low mean Leo among AGNs 
detected by PACS at 1.5 < z < 2.5. 

PACS-undetected AGNs and inactive galaxies also have 
comparable mean FIR luminosities. From the detection rates 
comparison in the inset panels, the general detection rate of 
AGNs is significantly higher than that of inactive galaxies among 
both isolated and interacting systems. However, the difference 
in detection rates between AGNs and inactive galaxies is not 
clearly dependent on the IM of the hosts, suggesting that star- 
bursts and nuclear activity are not any more coevally connected 
in mergers than in isolated galaxies. 

4.3. SFR as a function of Clumpiness Metric 

We have shown in Section 3.5 that AGNs are more likely to be 
in clumpy galaxies, at least at z < 1.5. In Figure 15, we ex- 
plore trends between FIR luminosity and the visual dumpiness 
of galaxies. As before, we look first at the mean Leo of AGNs 
and the inactive control sample in bins of redshift and CM (top 
row of the Figure). Firstly, there is a trend towards mildly ele- 
vated Leo among clumpy galaxies at all redshifts, which is ex- 
pected since clumps are a signature of unstable star-forming 
disks (Dekel et al. 2009; Ceverino et al. 2010). At z < 1.5, the 
enhancement in mean Leo for AGNs is small or non-existent 
among clumpy systems, but appears to be larger for smooth sys- 
tems. The differences between the mean Leo of AGNs and in- 
active galaxies as a function of CM reduces and possibly goes 
away by z ~ 2. 

The two lower rows of Figure 15 compare Leo distributions 
of the two populations. To allow sufficient number statistics, we 
use coarser bins in CM than in the upper panel: 0.0 <IM< 0.3 
indicates smooth galaxies, while 0.3 <IM< 1.0 includes both 
mildly and strongly clumped galaxies. AGNs and inactive galax- 
ies show fairly similar distributions in Leo in both CM bins, both 
among PACS-detected and -undetected systems. In the high red- 
shift bin, the AGNs are mildly weaker in mean star-forming 
properties compared to the inactive population. A comparison 
of PACS detection fractions shows again that AGNs tend to be 
detected a little more often than the control sample. There is a 
mild indication that the differences in detection rates are smaller 
among clumpy systems than among smooth systems. 

5. Discussion 

5. 1. Differences in structure and star-formation between 
AGNs and inactive galaxies 

In earlier sections, we analyzed structural and SFR patterns of X- 
ray selected AGN host galaxies and compared them to equally 
massive inactive galaxies. Broadly, AGNs are structurally sim- 
ilar to non- AGNs, whether in terms of galaxy light profiles, 
dumpiness or the incidence of interactions/mergers. Both sets 
of galaxies are primarily isolated, smooth and moderately bulgy 
systems. AGNs of all structural categories are more likely to 
detected in the FIR, although their SFR distributions, tracked 
using Leo distributions and stacks, are not vastly different from 
inactive galaxies (a more extensive treatment can be found in 
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Rosario et al. 2013b). Despite this, we do uneover some signifi- 
eant differenees through the use of a eareful eomparison of sta- 
tistieal distributions. 

First, AGNs show little systematie variation in their rest- 
frame optieal light profiles despite the obvious ehanges in the 
profiles of non- AGNs with redshift. At z ~ 1, massive inaetive 
galaxies are primarily bulgy, with a typieal Sersie index of 3, 
similar to the AGNs. By z ~ 2, inaetive galaxies beeome mostly 
disk-like and show Sersie indiees peaking at ^ = 1 . In eontrast, 
the typieal Sersie index of AGNs ehanges little between z = 2.5 
and z = 0.5. Comparisons of axis ratio distributions also high- 
light the rounder profiles of AGN hosts at z ~ 2. 

Comparing H-band and z-band light profile distributions re- 
veals that these differenees are produeed by a eentral light ex- 
eess in AGNs that is at least as red as the rest of the galaxy. 
Rest- frame UV light profiles of AGNs and non- AGNs at these 
redshifts are similar to their rest- frame optieal light profiles. 
This suggests that the star-forming disk eomponent, whieh likely 
dominates the rest-frame UV light, is eomparable in shape in 
both populations of galaxies, whereas the eentral red exeess is 
more pronouneed in AGNs. At z > 1, the higher Sersie AGN 
hosts are also more luminous, pointing to a possible relation- 
ship between the origin of the light exeess and the nuelear lu- 
minosity. While this exeess eould arise from reddened nuelear 
emission from an obseured aetive nueleus, the low luminosity of 
most X-ray seleeted AGNs in the CDFs implies instead that a 
more prominent stellar bulge may exist in AGN hosts (Seetion 
3.3.1 and Rosario et al. 2013a). Noting the important eaveats in 
this eonelusion outlined in Seetion 3.3.1, we proeeed with the 
eautious implieation of our results that the relationship between 
bulges and the existenee of super-massive blaek holes was al- 
ready in plaee at z ~ 2, although we refrain from any speeula- 
tion about the nature of that relationship or its evolution based 
on our present analysis. In the diseussion on AGN fueling modes 
(Seetion 5.3), we will eonsider the implieations of both alterna- 
tives: AGNs are in more bulge-dominant hosts at z ~ 2, or AGNs 
are in strueturally similar hosts to inaetive massive galaxies at 
z ~ 2. 

While the light profiles of AGNs are similar to inaetive 
galaxies at z ~ 1 (and probably not affeeted by AGN eontam- 
ination), their SFRs show different trends with Sersie index. In 
partieular, AGNs do not share the strong drop of SFR with n 
eharaeteristie of inaetive galaxies. An examination of the SFR 
distributions reveal that these high Sersie AGN hosts show both 
higher SFRs and a mueh higher FIR deteetion rate. This suggests 
that most bulge-dominated AGN hosts are indeed forming stars, 
while most bulge-dominated inaetive galaxies are fairly quies- 
eent. At z ~ 2, high Sersie AGNs still maintain high deteetion 
rates over inaetive galaxies, but the differenees in the mean SFRs 
are less pronouneed, probably beeause all massive galaxies, re- 
gardless of strueture, show larger SFRs at these redshifts. 

Using the visual elumpiness metrie (CM), we show that the 
majority of massive galaxies and AGN hosts at all redshifts are 
in relatively smooth galaxies. Despite this, we find a signifieant 
preferenee for AGNs to lie in elumpy hosts at z < 1. This pref- 
erenee weakens towards z ~ 2, while the relative fraetion of 
elumpy galaxies, both aetive and inaetive, inereases steadily. If 
the presenee of large elumps is related to the level of turbulenee 
or instability in galaetie gaseous disks, then the AGNs and inae- 
tive galaxies at z ~ 2 show the same level of sueh disturbanee, 
while at z < 1, AGNs appear to prefer elumpier galaxies and are 
more likely to be found in turbulent hosts. 

Aeross redshift, both aetive and inaetive elumpy galaxies 
have mutually eomparable SFRs. In addition, at low redshifts 


AGNs in smooth hosts display signiheantly enhaneed SFRs. 
Note that these smooth hosts inelude most of the high Sersie 
AGNs as well. Therefore, the large enhaneement in SFR ob- 
served among both bulgy and smooth AGN hosts are two sides 
of the same eoin, traeing mueh of the same population of host 
galaxies. 

Additionally, our study allows us to eompare the relative in- 
eidenee of interaeting or merging systems in AGNs and equally 
massive inaetive galaxies, using the visual interaetion metrie 
(IM) or the “multi-feature” merger subset. At z < 1, both trae- 
ers suggest a signifieant enhaneement of mergers among AGN 
hosts. The fraetion of interaeting systems among AGNs in the 
0.5 < z < 1.0 redshift bin is quite eonsistent with the frae- 
tions found for loeal hard X-ray seleeted Swift/BAT AGNs (18% 
- 25%; Koss et al. 2010). Towards z ~ 2, AGNs and inaetive 
galaxies have essentially identieal IM distributions. Therefore, 
the enhaneement in AGN aetivity in mergers seen in loeal and 
low redshift appears to beeome less pronouneed at higher red- 
shifts. 

Aeross IM, AGNs at z ~ 1 typieally show substantially 
higher mean SFRs than the eontrol sample. These differenees 
are driven by a eombination of mildly higher individual SFRs 
and a eonsistently higher FIR deteetion rate. The enhaneement 
in AGNs is eonsistent with the notion that nuelear aetivity is 
generally assoeiated with gas-rieh galaxies (Santini et al. 2012; 
Rosario et al. 2013b). We do not diseem any elear indieation that 
the SFRs of interaeting galaxies hosting AGNs are additionally 
boosted over isolated galaxies hosting AGNs, although this is 
predieted by most models of eo-eval SMBH fueling and star- 
bursts in galaxy interaetions. If there is sueh a seeondary en- 
haneement, it is minor when eonsidering the ensemble of vi- 
sually identified interaeting systems. On the other hand, this 
ensemble ineludes many gas-poor mergers, minor mergers or 
fly-bys, whieh do not neeessarily eonform to the predietions of 
merger simulations. A eareful treatment of merger samples, with 
additional elassifieation based on the gas eontent, mass ratios of 
the merging eomponents, as well as the stage of the interaetion, 
may reveal finer relationships. This may be possible in future 
work that uses the entire CANDELS area, but is beyond the sta- 
tistieal eapability of our analysis. 

5.2. Comparison with existing studies of X-ray selected 
AGNs 

A number of studies using various methods have explored the 
struetural properties of AGN hosts at intermediate and high red- 
shifts, sometimes using samples that overlap signiheantly with 
those used here. One advantage of our study is a eombination 
of visual and analytie measures of galaxy strueture over a wide 
range in redshift, whieh enables us to eompare our results to a 
broader subset of relevant work from the literature. 

We are generally eonsistent with earlier studies whieh used 
analytie measures of AGN host strueture, typieally at z < 1 
(Grogin et al. 2005; Pieree et al. 2007). Some of the differenees 
between the AGNs and inaetive galaxies in these earlier studies, 
sueh as higher early-type fraetions or bluer eolors among AGNs, 
are due to the improper ehoiee of eontrol samples, for example 
through matehing by a blue optieal luminosity rather than stel- 
lar mass. Sueh eontrol samples eontain a larger number of lower 
mass blue star-forming galaxies than ean be found among AGN 
hosts, whieh biases both struetural and SF measures. The need 
for a stellar mass-matehed eomparison sample is eritieal for a 
fair assessment of the AGN hosts (e.g., Villforth et al. 2014). 
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Kocevski et al. (2012) compared AGNs and mass-matched 
control galaxies at 1.5 < z < 2.5 in CANDELS/CDF-S using a 
visual classification scheme archetypical to the one in this paper. 
While we do not employ the visual separation into spheroids and 
disks as used in Kocevski et al. (2012), we have verified, using 
our visual classification catalogs, that our results are completely 
consistent with this earlier work (Appendix B). Among AGNs 
at these redshifts, most host galaxies are classified as disks, with 
a low merger/disturbed fraction and only a minor enhancement 
in the visual spheroidicity over inactive galaxies. At first glance, 
this seems to be at odds with the high Sersic indices seen among 
AGNs at z ~ 2. These differences stem primarily from the higher 
sensitivity of visual classifiers to the appearance of disks and 
a rather weak ability to discern variations in steepness of light 
profiles among spheroidal systems. An elongated n = 2.5 galaxy 
will be classified as a disk as easily as an elongated n = \ galaxy. 
Only GALFIT light profile modeling can adequately reveal the 
central light excess we find in AGN hosts at z ~ 2. A similar 
limitation will apply to visual studies of lower redshift AGNs 
from optical images such as Cistemas et al. (201 1), although the 
higher resolution of the HST/ACS imaging could potentially re- 
veal bulges more effectively. Visual methods are considerably 
more sensitive to distortions, asymmetries and disturbances in 
galaxies. 

From HST/WFC3 early release science imaging in GOODS- 
S, Schawinski et al. (2011) determined that AGN hosts at 1.5 < 
z < 3 were mostly in disky hosts, with typical GALFIT-based 
Sersic indices peaking at n = 1 . This result is at odds with our 
finding that AGNs at 1.5 <z< 2.5 have typical w of 2.5. One ma- 
jor difference lies in the fact that the AGNs in Schawinski et al. 
(2011) were fit using a two component galaxy model, with the 
assumption that the central light excess in these objects arises 
due to nuclear contamination. We only fit our galaxies with 
a single model profile and our two band studies demonstrated 
that the central excess is red and likely due to a bulge. By 
subtracting away a large part of the central bulge component, 
Schawinski et al. (201 1) may be lowering the effective Sersic in- 
dices in their AGNs (see Appendix A for a demonstration of this 
effect). To better understand the differences between these re- 
sults, careful multi-component modeling of the WFC3 images of 
AGNs and inactive galaxies is needed. Simulations of the influ- 
ence of multi-component fits on the recovery of bulge and disk 
parameters in realistic galaxies is needed, akin to the approach 
taken by Simmons & Urry (2008) and Gabor et al. (2009) for 
AGNs imaged with ACS. 

On the topic of merger incidence, the higher merger rates 
we And among AGNs at z < 1 is consistent with the 

level of enhancement of AGN activity in low redshift merg- 
ers (Silverman et al. 2011; Sabater et al. 2013). In addition, the 
weakening of a merger connection at z ~ 2 is consistent with 
Kocevski et al. (2012). Our results at z < 1 are in some ten- 
sion with the HST/ACS study of XMM-COSMOS AGNs from 
Cistemas et al. (2011), which found no enhancement in the in- 
teraction fraction over the control sample of inactive galaxies 
in an overlapping redshift range, despite a similar statistical 
power to our own study. This may arise due to the higher me- 
dian luminosities of the XMM-COSMOS sample, roughly an 
order of magnitude greater than the CDFs at these redshifts. 
Alternatively, and more likely, the differences may stem from 
the different approaches used to constmct a control sample; 
we match galaxies by stellar mass, while Cistemas et al. (201 1) 
match in redshift and apparent F814W (I-band) magnitude, cor- 
recting for any emission from a nuclear point source. The use of 
a blue optical rest- frame band to match galaxies will allow lower 


mass star-forming galaxies with a low M/L to enter the control, 
leading to a different stellar mass distribution between AGNs 
and the control sample (Xue et al. 2011; Rosario et al. 2013a). 
While it is not immediately apparent how this can affect merger 
incidence, a careful assessment of matching criteria is required 
before these differences are to be understood. 

5.3. Insights into AGN triggering scenarios 

A radiatively luminous AGN, including one bright in the X- 
rays, is triggered by dense gas falling onto an accretion disk 
around an SMBH. A major area of inquiry in the held of 
active galaxies pertains to how gas gets from scales of the 
host galaxy down to the black hole. Violent processes, such 
as galaxy mergers, are very effective at stripping angular mo- 
mentum from gas, sending it into the centres of galaxies to 
produce dense compact stmctures which can fuel synchronized 
starbursts and luminous AGNs. However, secular processes can 
also bring gas into the vicinity of the SMBH, as evinced by 
the gas-rich circumnuclear environments of settled disk galax- 
ies such as our own Milky Way (e.g., Morris & Serabyn 1996; 
Kruijssen & Longmore 2013). This gas can amply fuel low and 
moderate luminosity AGNs (Hopkins & Hernquist 2006), and, 
depending on the mechanisms relevant for the small-scale in- 
flow of gas around the nucleus, could even fuel luminous phases 
such as quasi-stellar objects (e.g. Gabor & Boumaud 2013). 

The morphology of a galaxy is an aggregate of its evolu- 
tion over many epochs. Once a galaxy’s stars are redistributed 
by violent processes into a spheroid, this marker of its merg- 
ing and inflow history is preserved, even if further inflow of gas 
and stars may settle into a later-forming disk. Folding together 
the morphological information of an AGN host, as a measure of 
its integrated evolution, and its SFR, as a measure of its current 
evolutionary state, can potentially constrain the importance of 
various AGN triggering models. 

We use the results of our study at z < 1 to demonstrate, 
through a simple heuristic example, the way structure and SFR 
may be jointly employed in testing fueling mechanisms. More 
quantitative and discriminatory tests will involve a comparison 
to observable predictions from cosmological semi-analytic mod- 
els or hydrodynamic simulations. 

AGNs at all redshifts are found in all forms of massive 
galaxies: those with disks, bulges and pure spheroids; isolated, 
clumpy, interacting or merging. Clearly, neither purely violent 
processes or purely secular inflow can account for all AGN trig- 
gering. The median Sersic index of ~ 2.5 implies that the typical 
AGN host is a disk galaxy with a substantial bulge, as is also 
conflrmed by visual and analytic estimates from earlier studies 
(e.g. Grogin et al. 2005; Pierce et al. 2007; Gabor et al. 2009). 
The prominence of a spheroid or bulge is a signature that violent 
processes likely played a role in the structural evolution of these 
galaxies. But the presence of a cold disk component, as well as 
the similarity of the Sersic index distributions of active and in- 
active galaxies, implies that either those violent events directly 
fuel only a small fraction of AGNs, or that such violent mech- 
anisms do not preclude the simultaneous formation or preser- 
vation of a galaxy disk (Robertson et al. 2006; Hopkins et al. 
2009; Dekel et al. 2009; Ceverino et al. 2010). The typical SFR 
of AGN hosts in disk galaxies is comparable to those of in- 
active galaxies and consistent with that of the SF Sequence 
(Mullaney et al. 2012; Rosario et al. 2013b): AGN activity must 
be driven in these galaxies by secular processes that do not 
strongly disturb the star-forming equilibrium of their gas disks. 
Therefore, any consideration for AGN triggering scenarios must 
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be mixed, with both seeular and violent eomponents. Theoretieal 
insight suggests that this mix evolves with redshift in the low and 
moderate luminosity AGN population (Hopkins et al. 2013b). 

We eonsider a seenario where AGN triggering is not direetly 
linked to proeesses that govern the larger seale galaxy. In this 
view, the role of the outer galaxy and its environment is to sim- 
ply supply gas to the inner kpe around the SMBH, whieh may 
arrive either through violent torques and relaxation, or through 
longer seeular means. Onee there, small-seale physies will gov- 
ern the final infall of this gas to the aeeretion disk, modulating 
the duty eyele of the aetive phase. SF indieates the presenee of 
eold gas in a galaxy, whieh is why we find AGN aetivity is en- 
haneed in star-forming galaxies aeross almost all galaxy mor- 
phologies and forms. In galaxies with very low gas reservoirs, 
sueh as massive elliptieals, the ehanee of hosting enough gas in 
the eireumnuelear regions is quite low. Indeed, only the small 
subpopulation of early type galaxies with signifieant gas eon- 
tent eontain an X-ray bright AGN and also display deteetable 
star-formation. This explains the preponderanee of SF in early- 
type AGN hosts. Among disk-dominated galaxies, the presenee 
of elumps are a signature of higher gas fraetions and more tur- 
bulent disks, whieh is why AGN hosts are mildly more elumpy 
at z < 1 . 

However, this simple seenario does not adequately explain 
the higher ineidenee of mergers among AGNs at these redshifts. 
Therefore, a ehannel must exist whieh involves a direet eonnee- 
tion between merger-driven torques on galaxy seales and gas in- 
fiow on to the SMBH on nuelear seales (e.g. Hopkins & Quataert 
2010). An potential signature of this direet “violent” ehannel is 
an average enhaneement in the SFR of AGN-ho sting mergers, 
whieh has been reported in detailed studies of loeal galaxy pairs 
and post-merger remnants (Liu et al. 2011; Ellison et al. 2011, 
2013). Unfortunately, the number of inaetive visually elassified 
mergers in our sample is too small to obtain a strong FIR detee- 
tion at z < 1 , and we eannot test this notion adequately in this 
work, but will be able to address this more eompletely in future 
work using all five CANDELS fields. 

In the diseussion of the violent ehannel, the enhaneed SF in 
high Sersie AGN hosts eould instead be taken as evidenee for a 
starburst in these galaxies followed by rapid quenehing within 
the last 100 Myr, the eharaeteristie timeseale over whieh a bolo- 
metrie traeer sueh as the FIR retains memory of a burst of SF 
(Hayward et al. 2014). Studies of the star-formation history of 
loeal emission-line seleeted AGNs in spheroidal hosts indieate 
that they underwent a reeent starburst (Sehawinski et al. 2010), 
possibly related to the event responsible for their morphologieal 
transformation. However, a seenario where gas brought in from 
external aeeretion onto already quenehed spheroidal galaxies, 
whieh then inspires both a starburst and AGN aetivity, ean also 
explain these results (Simoes Lopes et al. 2007; Martini et al. 
2013), so this is not a strong eonstraint on the extent of this pro- 
eess. 

Nevertheless, if we assume that all elumpy, interaeting and 
high Sersie AGNs at 0.5 < z < 1.0 are triggered by a violent 
ehannel, we ean plaee a rough upper limit on the fraetion of 
AGNs fuelled this way by adding the fraetions of all three eate- 
gories of hosts. Clumpy hosts (CM> 0.5) and interaeting hosts 
(IM> 0.5) eombined (ineluding overlaps) aeeount for ^ 30% 
of the population, while ineluding post-merger elliptieal hosts 
(with w > 3.5) brings the fraetion up to ~ 60%. Thus, potentially, 
more than half of the population of low and moderate luminosity 
AGNs at these redshifts may be fueled by gas brought to their 
eentres by violent meehanisms. A K-S test indieates that these 
AGNs do not have X-ray luminosities that are different from the 


rest, implying, to the degree we ean test with our small sample 
size, potentially violently triggered SMBHs have aeeretion rates 
that are similar to those fueled by other proeesses. 

At higher redshifts (z ~ 2), AGNs and inaetive galaxies 
both show identieal fraetions of elumpy or interaeting/merging 
systems, while both AGNs and inaetive galaxies have simi- 
lar SFRs irrespeetive of strueture. In partieular, merging AGN 
hosts have identieal mean SFRs and FIR deteetion rates as the 
mass-matehed eontrol. These results point to a lesser role for 
the violent ehannel in direetly fueling AGN aetivity, sinee syn- 
ehronisation between star-formation and nuelear aetivity appears 
to be weak. Rather, the signifieantly higher gas surfaee densi- 
ties in high redshift galaxies (Taeeoni et al. 2010; Daddietal. 
2010a; Taeeoni et al. 2013), eoupled with a faster seeular evo- 
lution timeseale (Genzel et al. 2008), ean effeetively diseonneet 
galaxy seale evolution from nuelear fueling proeesses. If AGN 
fueling at high redshifts is primarily modulated by small seale 
proeesses, then how do we understand the higher Sersie, rounder 
light profiles for AGNs in sueh high redshift hosts? If this is 
primarily due to widespread eontamination of these profiles 
by reddened nuelear light, and the true host struetures are in- 
deed disk-like and eonsistent with other massive galaxies, then 
a strong ease ean be made that seeular fueling is the primary 
mode at z ~ 2 (Sehawinski et al. 2011; Koeevski et al. 2012; 
Rosario et al. 2012; Hopkins et al. 2013b). However, if there is a 
preferenee for AGNs to be found in bulgy galaxies (as energetie 
arguments seem to support), then the violent mode responsible 
for the formation of sueh bulges may still retain a eritieal role 
in fueling higher redshift AGNs. Further elarifieation will eome 
through the eareful assessment of the true light profiles of these 
AGN hosts in up-eoming CANDELS studies. 
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Appendix A: The inclusion of a central point source 
in GALFIT fits 

Nuelear aetivity ean systematieally alter the appearanee of AGN host galaxies 
by adding excess emission in the centre, which concentrates the light profile, in- 
creases the resultant best-fit Sersic index and makes the system appear more cir- 
cular. Even at the highest attainable resolutions, the nucleus is completely unre- 
solved in distant galaxies, appearing as a point source. It is reasonable, therefore, 
to model AGN host galaxies as a combination of galaxy stmctural components 
and a central PSF. However, bulges, which also make light profiles more cen- 
tral concentrated, cannot be easily distinguished from low levels of nuclear point 
source emission. Galaxy bulges at z ~ 2 have characteristic half-light diameters 
of ~ 2 kpc or 0.24” (e.g. Bruce et al. 2012), only slightly larger than the FWHM 
of the WFC3/F160W PSF (0.15”). Therefore, accurately distinguishing between 
bulges and point sources requires careful modelling of the instrumental PSF, and 
even then may still be inaccurate, since galaxies do not usually have regular light 
profiles. 

To test the performance of GALFIT fits with a PSF component for real z ~ 2 
galaxies, we compared single Sersic and two component (Sersic + central PSF) 
fits to the F160W images of galaxies in the CANDELS GOODS-S field. This 
field has the deepest NIR imaging among the CANDELS fields. For this exer- 
cise, we used a subset of galaxies from the full sample introduced in Section 2.4, 
which have been carefully fit with different light profile models for a study of the 
bulge properties of distant galaxies by (Lang et al. 2014). Details of the method- 
ology and fitting setup are published in that work. Only galaxies at 1.5 < z < 2.5 
were considered in this exercise, since it is at these redshifts where the differ- 
ence in light profiles between AGNs and inactive galaxies is most pronounced 
(Section 3.3). The galaxies were all selected to have with M^> 10^® M©. 

Each galaxy was first fit with a single elliptical Sersic profile, iterated nu- 
merous times over a grid of initial values to prevent the best fit from falling into a 
local minimum. The best fitting single Sersic profile was then used to initialise a 
second fit with the addition of a PSF component to the galaxy model. The centre 
of the PSF was restricted to within 2 pixels (0. 12”) of the centre of the Sersic pro- 
file, and its flux was initialised to 1% of the integrated magnitude of the galaxy. 
All galaxies, whether identified as inactive or active, were fit identically in this 
fashion. 

In Figure A. I, we compare the resultant Sersic index distributions of the 
X-ray AGNs to those of mass-matched inactive galaxies, where the distributions 
for the inactive galaxies have been determined using the bootstrapping procedure 
discussed in Section 2.4.1. In the left panel, we show the results from single 
Sersic fits. These are qualitatively similar to those shown in Figures 5 and 7, 
in that AGNs typically show significantly higher Sersic indices than inactive 
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galaxies. In detail, we find a slightly higher fraction of « = 4 AGNs than in the 
full sample. However, given the small number of AGNs in this subsample, the 
difference could arise from Poissonian variation. 

In the right panel, we plot the distributions of the best-fit Sersic index of 
the galaxy component in the two component fits. As expected, including a PSF 
component has lowered the resultant n of the AGN hosts, greatly increasing the 
fraction of disk-dominated systems. On the other hand, an comparison of the dis- 
tributions of the inactive galaxies between both panels in the Figure also demon- 
strates a reduction in the typical n for these galaxies as well. Since the inactive 
population is not expected to show widespread nuclear point source emission, 
we conclude that simple two component fits as used in this exercise also tends to 
remove light from potential bulge components, systematically leading to best-fit 
Sersic index distributions that are too disky. 



^ogHser - nriogHp,, (AB) 

Fig.A.2. Distributions of the flux ratio of the Sersie eompo- 
nent to the PSF eomponent (in magnitudes) from two eomponent 
GALFIT fits to AGNs and inaetive galaxies in the WFC3/F160W 
band. 1000 draws of a mass-matehed eontrol sample are anal- 
ysed to determine the Icrjlcr seatter in the n distributions for in- 
aetive galaxies, shown as dark/light grey zones in the histograms. 
Red open histograms show the distributions for the AGNs. The 
dashed blaek line is the distribution for the full inaetive galaxy 
sample for whieh two eomponent fits were performed, ineluding 
many lower mass galaxies not widely found among AGN hosts. 
This histogram has been sealed down to overlap with the peak of 
the AGN histogram, to allow a simple visual eomparison of the 
distributions. 

Another valuable test of these fits is a comparison of the fraction of light 
in the PSF and galaxy components for AGNs and mass-matched inactive galax- 
ies. The nuclear luminosities of X-ray AGNs will be higher than any weak or 
heavily obscured nuclear emission that may remain undetected in the inactive 
population. If the central excesses are indeed due to widespread nuclear point 
source contamination, we expect to find higher PSF fractions among AGNs. In 
Figure A. 2, we plot histograms of the difference between the H-band magni- 
tudes of the best-fit Sersic component and the best-fit PSF components from our 
fits. In addition to the AGNs (red) and mass-matched inactive galaxies (grey re- 
gions), we also show the full distribution for galaxies with two-component fits 
(black dashed line), which includes many more low mass galaxies than generally 
found among AGN hosts. The latter histogram has been scaled down in number 
to allow a visual comparison to the other distributions in the Figure. 

Despite their higher nuclear luminosities, the two component fits yield es- 
sentially indistinguishable PSF fractions in the AGNs and equally massive inac- 
tive galaxies. This is not simply due to limitations of the fits or local minima, 
since the overall distribution of PSF fractions includes a long tail to very low 
values not found among the more massive AGN hosts (or other massive inac- 
tive galaxies). These two component fits were initialised with 5 mag between 
the PSF and the Sersic component, but the best-fit difference is about 2.5 mag 
(10%). Considering that the excess central light is found to a similar degree both 
in AGNs and inactive galaxies, one may conclude that the central excess is likely 
to arise in a bulge rather than in a nuclear point source. This is consistent with ar- 
guments based on central colours and energetics from Section 3.3.1. We refrain 


from commenting on bulge fractions here - this requires detailed bulge-i-disk de- 
composition fits to both AGNs and inactive galaxies. 

Appendix B: Comparison with the CANDELS study 
of Kocevski et ai. (2012) 

Kocevski et al. (2012) studied structural differences between AGN and mass- 
matched inactive control galaxies at 1.5 < z < 2.5 in the CDF-S using essen- 
tially the same visual classification scheme as in this work. However, there are 
important differences related to sample selections that should be borne in mind 
when comparing our results to theirs. 

The X-ray source catalogs used in the two CANDELS studies are based 
on distinct reductions and source detection algorithms, leading to different sized 
parent samples. We employ the CDF-S 4Msec catalog of Xue et al. (2011) which 
consists of 740 sources, while Kocevski et al. (2012) use a more conservative 
catalog of 569 sources. The differences between these catalogs are discussed 
in Rangel etal. (2014). In addition, we adopt AGN-specific photometric red- 
shifts from Luo et al. (2010), while Kocevski et al. (2012) take redshifts from 
Wuyts et al. (2008) which are not optimised for AGNs. Despite these differ- 
ences, the number of X-ray AGNs at 1.5 < z < 2.5, after the application of a 
lower Lx limit, is nearly the same in both studies (^ 70). Unlike us, however, 
Kocevski et al. (2012) do not restrict their sample to the GOODS-MUSIC foot- 
print or apply cuts in mn and M^. Our sample of AGNs for the visual classifica- 
tion analysis at 1.5 < z < 2.5 is 55, in total after cuts, a reduction of 25%. This 
smaller sample shares the same uniform analysis, quality checks and photomet- 
ric selections of the full GOODS-MUSIC dataset, enabling a consistent analysis 
of stellar mass and other galaxy properties between AGNs and inactive galaxies 
in this work. 

We can compare the fractions of disks and spheroids in our sample of AGNs 
with those reported by Kocevski etal. (2012). From their Table 1, ^ 80% of 
AGNs are classified to have a visible disk or spheroid. In contrast, 100% of our 
AGNs have one of these components. This is due to the mn < 24.5 cut applied to 
the galaxies in our visual classification catalog; we have less AGNs in our sam- 
ple, but all have accurate structural assessments. As a consequence, our fractions 
of disks and spheroids will necessarily be higher than those in Kocevski et al. 
(2012), simply due to our different sample sizes. Therefore, we scale our frac- 
tions down by a factor of 1 .25 to ease the comparison. 

In our AGN sample, we find disk galaxy fractions of 55^^% of which 
13+5% are pure disks, with no reported spheroid component. These fractions 
maybe compared to 51+^ and 17+4% respectively from Table 1 of Kocevski et al. 
(2012). We also find pure spheroid fractions of 25+5 % compared to 26+^ from 
Kocevski et al. (2012). As for morphologically disturbed systems (Section 3.4), 
the fractions of visual disks and spheroids are completely consistent in both 
CANDELS works. Our ability to reproduce the results of Kocevski et al. (2012), 
despite the differences of approach and numbers of classifiers, highlights the sta- 
bility of the CANDELS visual classification scheme. 

This being said, we prefer in this work to use an analytical measure of the 
galaxy light profile rather than visual measures of diskiness. As stated in Section 
5.2, a visual classifier has difficulty discriminating between subtle variations in 
the light profile gradient. For e.g., inactive galaxies at z ~ 2 that are classified 
visually as having dominant disk components exhibit a range in Sersic index 
of 0.7-3. 2 (80th percentile), while the range is 1.7-6. 6 for spheroid-dominated 
galaxies. 
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Fig. 13. Mean 60 /um monochromatic luminosities (Leo) and luminosity distributions of AGNs and inactive control galaxies in both 
GOODS fields combined, as a function of H-band Sersic index (n). Panels left to right span three distinct bins in redshift. In the 
top row, the mean Leo from combinations of detections and stacks are compared in bins of n. The X-axis error bars show the extent 
of these bins. The errors on the AGN measurements (red points) are determined from bootstrapping into the AGN subsample in 
each bin in redshift and n. The errors on the inactive galaxy measurements (black points) come from the analysis of an ensemble of 
mean Leo from 100 random draws of mass-matched control galaxies. In the lower two rows, each panel is split into two subpanels. 
Leo distributions and PACS detection fractions for AGNs and inactive galaxies are plotted in the left subpanels for two coarse bins 
in Sersic index. Histograms show the Leo distributions of PACS-detected AGNs and inactive galaxies. The 100 random draws of 
mass-matched control galaxies are used to determine the lcr/2cr scatter in Leo for inactive galaxies, shown as dark/light grey zones in 
the histograms. The red histograms are the Leo distributions for AGNs. The mean Leo corresponding to these histograms are plotted 
for PACS-detected AGNs as red star points and for inactive galaxies as black star points. Error bars on the latter are the rms scatter 
of the mean of the control samples. The mean Leo of PACS-undetected AGNs (red circle points) and inactive galaxies (black circle 
points) are also compared similarly. PACS detection fractions are shown in the right subpanels. Red cross points denote AGNs and 
black cross points denote inactive galaxies, with rms scatter shown as vertical error bars on the latter points. See Section 4.1 for a 
discussion. 
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Fig. 14. Mean 60 //m monochromatic luminosities (Leo) and luminosity distributions of AGNs and inactive control galaxies in 
GOODS-S, as a function of visual Interaction Metric (IM). Panels left to right span three distinct bins in redshift. In the top row, the 
mean Leo from combinations of detections and stacks are compared in bins of IM. In the lower two rows, each panel is split into 
two subpanels. Leo distributions and PACS detection fractions for AGNs and inactive galaxies are plotted in the left subpanels for 
two coarse bins in IM. PACS detection fractions are shown in the right subpanels. Details of the plot are identical to those of Figure 
13, except for a difference in the structural measure. See Section 4.2 for a discussion. 
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Fig. 15. Mean 60 /um monochromatic luminosities (L6o) and luminosity distributions of AGNs and inactive control galaxies in 
GOODS-S, as a function of visual Clumpiness Metric (CM). Panels left to right span three distinct bins in redshift. In the top row, 
the mean L^o from combinations of detections and stacks are compared in bins of CM. In the lower two rows, each panel is split 
into two subpanels. L^o distributions and PACS detection fractions for AGNs and inactive galaxies are plotted in the left subpanels 
for two coarse bins in CM. PACS detection fractions are shown in the right subpanels. Details of the plot are identical to those of 
Figure 13, except for a difference in the structural measure. See Section 4.3 for a discussion. 
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Fig.A.l. Distributions of the Sersic index {n) from GALFIT fits to AGNs and inactive control galaxies in the WFC3/F160W band. 
1000 draws of a mass-matched control sample are analysed to determine the lcr/2cr scatter in the n distributions for inactive galaxies, 
shown as dark/light grey zones in the histograms. Red open histograms show the distributions for the AGNs. Left: Results from 
fits of a single Sersic elliptical model. Right: Results from fits with two components - a Sersic elliptical model and a central point 
source. The inclusion of a point source results in lower Sersic indices for both AGNs and inactive galaxies. 
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